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Abstract of the Dissertation 

Single Electron Transverse Momentum and 
Azimuthal Anisotropy Distributions: 
Charm Hadron Production at RHIC 

by 

Xin Dong 

Doctor of Philosophy in Physics 
University of Science and Technology of China, Hefei, 2005 

Quantum Chromodynamics (QCD) is a basic gauge field theory to describe strong 
interactions. Lattice QCD calculations predict a phase transition from hadronic matter to 
a deconfined, locally thermalized Quark-Gluon Plasma (QGP) state at high temperature 
and small baryon density. Plenty of exciting results from RHIC experiments in the first 
three years have demonstrated that a hot dense matter with strong collective motion which 
cannot be described with hadronic degrees of freedom was created at RHIC. Charm quarks 
are believed to be mostly created from initial gluon fusion in heavy ion collisions. Since 
they are massive, charm hadrons are proposed to be ideal probes to study the early stage 
dynamics in heavy ion collisions. 

We provide here an indirect measurement of charm semi-leptonic decay. Single electron 
transverse momentum (p^) distributions from 200 GeV d + Au , p + p coUisions and 62.4 
GeV Au + Au coUisions, and single electron azimuthal anisotropy (^2) from 62.4 GeV Au 
+ Au collisions are presented. 

Electron identification is performed with the combination of the prototype Time-Of- 
Flight detector (TOPr) and the ionization energy loss (dE/dx) in the Time Projection 
Chamber (TPC). Photonic background electrons are subtracted statistically by recon- 
structing the invariant mass of the tagged and every other partner candidate e^. 
Partner track finding efficiency is estimated from Monte Carlo simulations. More than 
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~ 95% of photonic background (photon conversion and 7r° Dahtz decay) can be sub- 
tracted through this method. The non-photonic electron px spectrum is extracted. In 
Au -I- Au coUisions at sJsnn — 62.4 GeV, due to low yield of charm quarks, only inclu- 
sive electron and photonic electron pt spectra are presented to illustrate the feasibility 
of this method in Au -|- Au collisions. Electron azimuthal anisotropy in Au -|- Au 62.4 
GeV is calculated from the event plane technique. Photonic background in each bin 
is reconstructed using the invariant mass method. V2 of inclusive electrons and photonic 
electrons is presented. 

The non-photonic electron pt spectrum is consistent with the direct open charm re- 
construction spectrum va. d -\- Au collisions. The total charm production cross section 
per nucleon-nucleon collision is calculated from the combined fit of the spectrum 
and non-photonic electron spectrum in o? -|- Au collisions. The result is da^^^/dy — 
0.30±0.04(stat.)±0.09(syst.) mb, which is significantly higher than the Next-to-Leading- 
Order perturbative QCD (pQCD) calculations extrapolated from low energy data points. 
The total charm cross section measurement is crucial to investigate J/ip production mech- 
anism in Au + Au collisions, which may be a robust signature of QGP formation. Not 
only the total yield is larger, but also the pt spectrum is harder than those from pQCD 
predictions, indicating a possible unusual charm fragmentation function at RHIG energy. 
In Au -I- Au 62.4 GeV, the photonic electron spectrum agrees with the inclusive electron 
spectrum well. The V2 of photonic electrons also agrees with that of inclusive electrons. 
These are consistent with the expected low charm yield at this energy. The success of 
this technique in Au -|- Au 62.4 GeV data makes us confident about extracting both the 
non- photonic pr spectrum and V2 in the coming large 200 GeV Au -|- Au data set. 

Many open charm measurements can be carried out with a full barrel TOF and heavy 
flavor tracker upgrade. These upgrades will allow us to study the details of open charm 
collective motion, eUiptic flow, correlations and so on with much better statistics. These 
two sub-detector upgrades will also allow us to reconstruct vector mesons (cj, (p etc.) from 
the low mass di-electron invariant mass distributions, which may offer us indications of 
possible chiral restoration issues during a possible phase transition at RHIG. 
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CHAPTER 1 



Introduction: Quantum Chromodynamics and Heavy 

Ion Collisions 

1.1 Quantum Chromodynamics 

From the lepton-nucleon Deep Inelastic Scattering (DIS) experiments in the late 1960's, 
it was formed that hadrons have structure and the quark parton model was exposed. The 
introduced sub-nucleon structures include quarks which are constituents of hadrons, and 
gluons which propagate interactions between partons. This kind of interaction, namely 
the strong interaction, is one of the four fundamental interactions in the nature. Quantum 
ChromoDynamics (QCD) has been established since the 1970's to describe strong inter- 
action, and together with the unification of electroweak theory, composes the Standard 
Model (SM), which is a reasonably successful model to describe all interactions except 
gravity. 

QCD |DKS03j is based on the gauge group SU{?>)c, with gauge bosons - color octet 
gluons for factors and a unique group coupling constant Qs- The subscript C denotes 
the quantum number - color, which is an exact symmetry. Quarks belong to a color 
triplet representation in this symmetry, but hadronic states are assumed to be color 
singlets in QCD. Owing to the non-abelian character of the color group, the invariant QCD 
Lagrangian^ requires gauge (gluon) self-interactions, which do not appear in Quantum 
ElectroDynamics (QED) - the gauge theory describing electromagnetic interaction. 

Generally speaking, QCD is a non-perturbative gauge theory in most cases. It can be 
^See Appendix A for QCD Lagrangian 
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calculated using a computer-assisted method - Lattice QCD |Gup98| . In this calculation, 
the spacetime is discretized and replaced by a lattice with lattice spacing equal to a. 
The quark fields are only defined at the elements of the lattice and the gauge fields are 
defined on the links of the lattice. The action is rewritten in such a way that the limit 
a — i> formally gives the original continuous action. Lattice QCD has been widely used 
for reliable QCD calculations. The precision of Lattice QCD calculations are limited by 
the lattice spacing or the computing power. 

1.1.1 QCD running coupling constant ag 

The renormalized QCD coupling shows a scale dependent coupling as{fi) (running 
coupling), similar to that in QED. However, the QED running coupling increases with 
energy scale, while the gluon self-interactions lead to a complete different behavior in 
QCD. as(/i) can be written as: 

"^(^) - 4^ - , , , t:.. , (1-1) 

When Po > 0, this solution illustrates the asymptotic freedom property: — *• as /i — 
oo, which means QCD can be calculated perturbatively in high momentum transfer or 
short distance approach. On the other hand, this solution also shows strong coupling 
at /i ~ Aqcd, so QCD is non-perturbative in this case, as needs to be determined 
from experiment. The world averaged as at the fixed-reference /iq = Mz is as{Mz) = 
0.1 187± 0.002 |Eid04j . and the QCD scale Aqcd ~ 200 MeV. Fig. ITT] shows the measured 
as at different momentum transfer scale /x compared with Lattice QCD calculations. 

1.1.2 Perturbative QCD (pQCD) 

At sufficiently high /i, where as is sufficiently small, physics quantities, such as cross 
sections, can be calculated to a truncated series, known as Leading Order (LO), Next- 
to-Leading Order (NLO) etc.. There are plenty of experiments on high energy processes 
which offer quantitative tests of pQCD. Due to the complexity of hadronic processes, the 
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Figure 1.1: Measured QCD running coupling constant a<j from different experiments 

compared with Lattice QCD calculations. 

test for pQCD has more difficulties than that for QED. 

Assuming factorization, the cross section of a process A + B ^ C + ... can be written 

as: 

(TAB^C = fa/A{Xa, fiF)fb/B{Xb, fij,) ® ^ab^c{s, ® D^^ciz, fl%) (1.2) 

Only the middle term aab^c can be calculated in pQCD from Feynman diagrams. The 
first term fa/A{xa, f^'p) or fb/si^b, fJ''F) is the hadron Parton Distribution Function (PDF) 
and the last term Dc->c{z, fi'p) is the Fragmentation Function (FF) that describes the 
transition from a parton to a hadron. For leptons, these two terms do not contribute in 
this formula. Hence, we can measure PDFs through lepton-nucleon DIS interactions and 
FFs through high energy e+e~ collisions, fi^i is the renormalization scale, originating from 
the need to regularize divergent momentum integrals in calculating high order diagram 
loops. yUi? is the factorization scale, at which the parton densities are evaluated, s is the 
partonic center of mass energy squared. From this formula, we can see it is complicated 
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to determine the expected hadron production cross section in hadron-hadron collisions. 

Heavy quark (c, b) production, due to large masses, is believed to match to pQCD 
prediction better than light quark production. And because they cannot be produced 
through the initial light hadron fragmentation, the FF part is irrelevant to the total 
production cross section of heavy quarks. Hence the measurement of the total heavy 
quark production cross section offers a powerful test of pQCD. 

Let us take the calculation of heavy quark pair QQ cross section in p + p collisions as 
an example |Vog02 . At LO, heavy quarks are created by gg fusion and qq annihilation. 



while at NLO, qg and qg scattering is included. To any order, the partonic cross section 
can be calculated as an expansion in the power of a^, assuming fiji = fip = fi: 

2 / \ oo k 2 

a,,{s,ml^^') = ^ 5:(4vra.(^))^ (1-3) 

A:=0 /=0 ^ 

where ?7 = s/Arn^ — 1 and fij'^\f]) is called the scaling function, obtained from the 
calculation of Feynman diagrams for each order. The total production cross section can 
be obtained from the above partonic cross section: 

app{s,ml)= [.—6{xiXj-T)fi/p{xi,fi^)fj/p{xj,fi'^)crij{T,ml,fi^) (1.4) 

i,j=q,q,9 s 

Numerical results of theoretical calculations show strong dependence on the truncation 
of the series, the selection of parameters fip, and niQ etc. |MNR93| Vog02 IRPOSj . 
Hence, theoretical predictions still have large uncertainties. 

Experimentally, measurement of charmed hadrons is difficult due to their short life- 
time (cr(Z)°) = 124 /im), low production rates, and large combinatoric background. 
Both direct reconstruction through the hadronic channel and indirect measurement of 
the semi-leptonic decay of charmed hadrons were preformed in the previous measure- 
ments. In the low center-of-mass energies ( ^ 40 GeV ), measurements were done on 
fixed targets jTav87[ IBar881 lAmmSSL IKodQlllXK^ . At ~ 52 - 63 GeV (ISR), the 
measurements were done from DK^ production and ee, e/x pair correlation, and the results 
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showed inconsistency between different publications |Tav87j . At higher energies, the UA2 
experiment measured single electron distributions at 630 GeV p + p collisions |Bot90j . 
However, due to large uncertainties of the decay from charmed hadron to electrons, these 
measurements included large errors. The CDF II collaboration made a direct measure- 
ment of open charm hadrons at ^/s = 1.96 TeV p + p collisions |Aco03j . but the spectrum 
only covers the high px region. Recent pQCD calculations seem to reproduce the spec- 
trum at CDF. No charm measurement in elementary p + p collisions at RHIC energy 
was done and theoretical predictions differ significantly at this energy due to the large 
uncertainties of many parameters. So charm production at the baseline elementary p + p 
collisions at RHIC energy is helpful. 

1.1.3 Confinement and chiral symmetry breaking 

Since quarks have color quanta, while hadrons are color-neutral to us, quarks must be 
confined within hadrons. This can also be explained from the QCD coupling as- When 
two quarks are separated to a large distance, which corresponds to a small energy scale, 
the coupling becomes strong, i.e. intuitively, more and more self-coupled gluons hold the 
quarks not to be isolated. It is quite different in QED because there is no self coupling 
between photons, so that we can observe isolated electric charges. 

In the absence of quark masses, the QCD Lagrangian can be split into two independent 
sectors: the left- and right-handed components |Pic95j . This Lagrangian is invariant under 
chiral symmetry transformations then. This symmetry, which is the extension of classical 
SU{3), is a global SUiinf) x SUR^Uf) symmetry for Uf massless quark flavors. However, 
it is spontaneously broken in the vacuum in the Nambu-Goldstone's way to realize this 
symmetry and this breaking gives rise to {nj — 1) massless Goldstone particles. Thus we 
can identify the vr, K, rj with the Goldstone modes of QCD: their small masses being 
generated by the quark-mass matrix which explicitly breaks the global chiral symmetry 
of the QCD Lagrangian^. 

^See Appendix A for chiral symmetry and effective lagrangian approach. 
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1.1.4 QCD Phase transition 

QCD matter is mostly observed as nuclei or hadron gas in our current condition. 
With sufficient temperature and energy density, QCD predicts a phase transition to a 
new matter, named Quark Gluon Plasma (QGP), with new (color) degrees of freedom 
(d.o.f.). In this new phase, quarks and gluons are liberated from hadrons, and can move 
around in a larger distance rather than confined in hadrons, which is called deconfinement. 
Meanwhile, the broken chiral symmetry in normal QCD matter will be restored and 
consequently, masses of scalar mesons and vector mesons will decrease |Kar02aj . Lattice 
QCD calculations provide quantitative predictions on this phase transition: the critical 
temperature of this phase transition is ~ 150 — 180 MeV, and the energy density at the 
critical point is ec{Tc) ~ 1 — 3 GeV/fm'^ (~ 0.17 GeV/fm^ for nuclear matter) |Kar02bj . 
The appearance of these color d.o.f. can be illustrated by a sharp increase in pressure 
with temperature, shown in Fig. 11.21 |KLP00j . 

1- 
Q. 
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12 3 4 

T/T, 

Figure 1.2: The evolution of p/T^ with the increase of temperature T for 3 different 
flavor configurations. The arrows indicate the SB limit for each case. The insert plot 
shows the ratio of p/psB with function of T. 

The arrows indicate the Stefan-Boltzmann limits, which are for the systems with 
massless, non-interacting quarks and gluons. The similarity of the three curves in the 
insert plot of Fig. 11.21 illustrates that besides the effect of quark masses, there should 
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be interactions in the newly formed system, which is different from the original QGP 
scenario (weakly interacting), but is already demonstrated by experimental results (see 
next section). 

Lattice QCD calculations of the potential between two heavy-mass quarks also offer 
evidence of deconfinement. 

Lb 2 
> 

1.5 
1 

0.5 


1 2 3 4 5 

Figure 1.3: Lattice calculations for the heavy-mass quark potential in different 
temperature cases. The band depicts the Cornell potential of V{r) = —a/r + ar with 

a = 0.25 ±0.05. 

Fig. 11.31 shows a recent calculation of the heavy-mass quark-antiquark pair Cornell 
potential in different temperature conditions |KLP00j : with the increase of temperature, 
the rampart of the potential between two quarks, which causes confinement, will bend 
down and thus liberate quarks from the trap. In addition, the continuous bending without 
sudden change indicates a crossover transition at high temperature and vanishing net 
quark density. 

1.2 Heavy Ion Collisions 

Experimentally, to search for this new kind of matter, a large amount of energy needs 
to be packed into a limited space volume. Heavy ion collisions have been proposed as a 
more effective way because the initial energy density increases as a power law function 
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with the atomic number while only logarithmically with collision energy |Lin96j . Since 
the 1970's, physicists from BEVALAC at LBL, SIS at GSI, AGS at BNL, SPS at CERN 
and RHIC at BNL etc |Sto04j have been trying to reach and cross the phase transition 
boundary in the laboratory through relativistic heavy ion collisions. 



Quark Gluon Plasma 

point 




0.0 0.2 0.4 0.6 0.8 1.0 



Baryonic Chemical Potential |ig (GeV) 

Figure 1.4: QCD phase diagram. The grey shadows depict first-order phase transition 
boundaries. The red dots depict the critical points and "3cr point" is calculated from 
Lattice QCD |FKn2j . The blue dots depict the positions of several colliders from 

statistical fit |Bra01j . 

Fig. 11.41 shows the QCD phase diagram and approaching transition boundary of the 
colliders. 

For the past four years the Relativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Lab (BNL) has conducted very successful runs. Plenty of exciting physics results 
reveal that the matter created at RHIC is quite different from what we observed before: 
It cannot be described by hadronic degrees of freedom and demonstrates many of the 
signatures from a QGP scenario. These measurements provide strong hints for the dis- 
covery of QGP |Adc04j . Some of the key measurements will be discussed in the following 
sections'^. 



•^Useful kinematic variables are defined in Appendix B. 
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1.2.1 Nucleon stopping power and initial energy density 



One of the most fundamental quantities we need to investigate is whether the initial 
deposited energy is enough to cross the energy density threshold for QGP formation. 
Because baryon number is conserved and the rapidity distributions are only slightly af- 
fected by rescattering in the late stage of collisions, the measured net baryon {B — B) 
distribution can reveal the energy loss of initial participants and allow us to estimate the 
degree of nucleon stopping power. Fig. II. 51 shows rapidity loss {5y) = {y) —yp for different 
energies from AGS to RHIC |Bea04j . Using the data point at RHIC, one can estimate 
that 73 ± 6tll GeV of the initial 100 GeV per participant is deposited and available for 
excitations. 



A 

to 

V 

CO o 

CO 3 
O 



P 2 
Q. 
CC 

cc 

1 



50 



-dN/dy Net Baryon A 




1 ' 

r 


1 ' 


— momgauss 1 \ 

— pol6 A'' 












n 


12 3 4 


5 




1 




X- 



A E917 
T E802/E866 
■ NA49 
• BRAHMS 



_L 



1 2 3 4 5 

Beam rapidity yp 

Figure 1.5: The rapidity loss as a function of beam rapidity. The shadow indicates an 
unphysical region and the dashed line depicts the phenomenological scaling 
{5y) = 0.58|/p. The insert plot shows the measured data points of net baryon 
distribution from BRAHMS and different extrapolations to full rapidity. 



The initial Bjorken energy density |Bjo83| can be calculated using 

1 dEr 



eBj 



A±T dy 



:i.5) 



where r is the formation time and A± is the nuclei transverse overlap region area. 
The PHENIX dEx/dy measurement |Adl04aj indicates an initial energy density of ~ 5 
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GeV/fm^ (r ^ 1 fm/c, A± = tcR"^, R ^ 1.2A^/^ fm) for central Au + Au collisions at 
RHIC, well above the expected critical energy density Cc ~ 1 GeV/fm^. 

Thus, the initial condition in Au + Au collisions at RHIC is believed to be capable of 
forming the new partonic matter - QGP. 



1.2.2 Jet quenching 

In heavy ion collisions, high {pr >~ 5 GeV/c) particles are believed to be pro- 
duced mainly from the initial QCD hard-scattering processes |Adin2aj . These energetic 
particles can be used as unique probes by studying their interactions with the medium. 
Experimentally, the nuclear modification factor, the difference between the spectrum in 
A+B collisions with respect to (w.r.t.) a. p + p collision reference, is extensively used. It 
is defined as: 

(PNABldprdy 

Rab{Pt) = 7^^-72 / , , 1-6 

where Tab = {^bin)/^^^^^ is the nucleus overlap function, calculated from a Glauber 
model |Won94j . Ni,in represents the number of binary collisions in a nucleus- nucleus 
collision. The experimental result |Ada03bj in Fig. 11.61 shows that there is a strong 
suppression relative to the binary scaling expectation at high p^ in the central Au + Au 
collisions — jet quenching. But this suppression is not seen m. d + Au collisions, the 
control experiment, which suggests the suppression in central Au + Au is due to the final 
state interactions rather than initial state effect and thus a very dense matter must be 
created in central Au + Au collisions at RHIC. 

This result has also been illustrated by a dihadron azimuthal angle correlation study. 
Fig. 11.71 shows the associated hadrons {pt > 2 GeV/c) azimuthal distribution relative to 
a triggered hadron {p^ > 4 GeV/c). The enhanced correlation at A0 ~ 0, which means 
the pair is from a single jet, was observed in p + p , d + Au and Au + Au collisions. The 
pair from the back-to-back jet correlation at A0 ~ vr only appears in p + p and d + Au 
while it almost completely disappears in central Au -|- Au collisions. 

Due to momentum conservation, the disappearance of away-side fast partons must 
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result in an excess of softer emerging hadrons. Fig. II. 81 shows the centrahty dependence of 
the {pt) of the associated away-side charged hadrons (threshold lowered to 0.15 GeV/c), 
compared with that of inclusive hadrons, illustrating the above point j Wan04j . This 
also offers a hint of the attainment of thermalization via the frequent soft parton-parton 
interactions in the early collision stages. But how strong those partons interact is still a 
crucial open question, that needs to be answered quantitatively to address the evidence 
of early thermalization of the system. 

1.2.3 Collective motion 

Hadron spectra are useful tools to study the properties of the bulk system and specified 
particles can be used to probe different stages after the heavy ion collisions. Fig. ll.9l shows 
the identified particle pt spectra measured in central Au+Au collisions |Ada04at IAda04b| 
IAdan4clin^^ . 
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Figure 1.9: Identified particle spectra in central Au+Au collisions at ^/snn = 200 GeV 
and the Blast Wave fit results. The BW fits were done for tt~ , K~, p simultaneously and 
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The plot shows the slopes of particle spectra changes for different particles (masses), 
indicating the strong collectivity of final state particles. The dashed lines depict the fit 
results from the Blast Wave thermal model |SSH93j . In thermal models, local thermal 
equilibrium is assumed and hence particles spectra only depend on the mass of particles 
and the temperature of system. "Blast wave" means that particles freeze out from the 
system surface simultaneously when dense matter becomes dilute enough. Under the 
assumption of simple cylindrical source and boost invariance in rapidity, there are only 
two parameters to describe the particles spectra: freeze-out temperature Tjo and aver- 
age transverse velocity (/3t)- Fig. II.IUI shows the fit results for different particles. For 
the simultaneous fit to tt^, K~ , p, stronger and stronger collectivity is observed from 
peripheral to central collisions. But the fit to Q and (p in central Au + Au collisions 
shows higher freeze-out temperature and lower transverse velocity, indicating those par- 
ticles leave the system at the earlier stage than stable hadrons. This is not surprising 
since those multi-strange baryons/mesons are expected to have much smaller hadronic 
scattering cross sections and thus the temperature from the fit to those particles may 
reflect the chemical freeze-out temperature Tch- This temperature is close to the critical 
temperature Tc, meaning the temperature of the system created in the collisions is greater 
than Tc and hence the phase transition may take place at RHIC energy. 

1.2.4 Transverse azimuthal anisotropy 

In non-central heavy ion collisions, the overlapping region of the two nuclei will form an 
anisotropy in coordinate space. Because of rescattering, the anisotropic pressure gradient 
will lead to an anisotropy in momentum space. The dynamic expansion of the system 
will wash out the coordinate-space-anisotropy, while the momentum-space-anisotropy will 
saturate during the evolution of the system |KH03j . The final state particle spectrum in 
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Figure 1.10: Blast-wave parameters Tfo vs. {(3t) contour plot from the simultaneous 
fits to stable hadrons (vr, K^p) spectra in Au + Au collisions and p + p collisions and 
separate fits to multi-strange hadrons VL{sss) spectra in central Au -|- Au 

collisions. The contours of vr, K,p fits in Au + Au are for peripheral collisions in the left 
and central collisions in the right. The contours of and VL are for central Au -|- Au 

collisions. 

momentum space can be expanded into a Fourier series as Eq. 11.71 



E 



dp^ 2'KpTdpTdy 



:i + ^2w„cos[n(0-^,p)]) 



n=l 



Vn = {C0s[n{(f) - "^rp)]) 



(1.7a) 
;i.7b) 



where ^^^p denotes the direction of the reaction plane. The first and second harmonic 
coefficients vi, V2 are called directed and elliptic flow. Due to the approximate elliptic 
shape of the overlapping region, the elliptic flow V2 is the largest harmonic observed in 
mid-rapidity. Because of the quenching of coordinate-space-anisotropy, elliptic flow can 
reveal early information about the system and because it depends on rescattering, elliptic 
flow is sensitive to the degree of thermalization of the system in the early stage. 

Fig. II .111 shows identified particle V2{pt) and the hydrodynamic model predictions for 
Pt < 2 GeV/c |Ada05j . In this low Pt region, V2 has larger values for lower mass particles. 
This mass ordering is reasonably described by the hydrodynamic models, which assume 
ideal relativistic fluid flow and negligible relaxation time compared to the time scale of the 
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compared with hydrodynamic model predictions. 

equilibrated system. The agreement implies early thermalization, i.e. strongly interacting 
matter with a very short mean free path dominates the early stages of the collisions. 
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Figure 1.12: Left: Identifled particle V2 up to intermediate pt ■ Right: After scaling 
the V2 and px with number of constituent quarks Ug, all particles fall onto one universal 

curve at px/^q > 0.6 GeV/c. 



Fig. 11.121 shows V2{pt) in a larger px range for different particles |Adl03a| lAda04bj . V2 
for all particles saturates above a certain px (~ 2 — 3 GeV/c). On the left panel, in addi- 
tion, particles show two groups on the plots: mesons and baryons. After scaling both V2 
and Pt with the Number of the Constituent Quarks (NCQ) in the corresponding hadron, 
all particles with px/ng > 0.6 GeV/c fall onto one universal curve except pions (due to 
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resonance decay effect |GK04| n3ES04j ). shown on the right panel. This meson/baryon 



grouping phenomenon was also observed in the nuclear modification factor Rcp at inter- 
mediate pt (1.5 < pT/(GeV/c)< 5) HHaOiEl EorOSl . Coalescence models pOMlMyna] 
which assume hadrons are formed through coalescing of constituent quarks provide a 
viable explanation for these observations. This indicates the flow developed during a 
sub-hadronic (partonic) epoch, and offers a strong evidence of deconfinement at RHIC. 

1.2.5 Heavy flavor production in HIC 

Heavy flavor hadrons are expected to be unique tools to probe the early stage in- 
formation in heavy ion collisions due to the following features. The creation of heavy 
flavor quarks is dominated by the initial gluon fusion processes and is negligible in the 
pre-equilibrium stage jEm96j. Systematic studies of charm production in p + p , and p 
+ A collisions have been proposed as a sensitive way to measure the PDF in the nu- 
cleon, and the nuclear shadowing effect |LG96j . Due to their heavy masses, energetic 
heavy flavor partons are expected to lose less energy than light quarks from gluon radia- 
tions when traversing the dense medium ("dead-cone" effect) |DK01j . At RHIC energies, 
due to the possible large production of charm quarks, coalescence processes might not 
be negligible, especially for closed charm production, indicating the standard J/ip sup- 
pression scenario might be invisible |(4Rnil ITSHniL l(4KLn21 lABROSj . Recent studies 
propose that charm quark flow not only offers a check of hydrodynamic models and/or 
coalescence models, but also could imply the thermalization of the light quarks in the 
system |DESn41 IBK(;n31 K;KRn4j . 

Theoretically, there are lots of calculations of heavy quark energy loss |DGW05l 
IDG04b| IDG03j , nuclear shadowing |LG96j , cronin enhancement |Vog03| , charmonium pro- 
duction |NL(;n8j and charm quark flow jLMOai IB(:;sn4j etc. for the collisions at RHIC. 
However, we have very few experimental results to test or prove these features so far. The 
first publication of charm production at RHIC energy was a single-electron measurement 
by the PHENIX collaboration at ^/s^n = 130 GeV in Au + Au collisions |Adc02j . They 
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claimed that charm production in heavy ion colhsions obeys the number of binary colh- 
sions scahng and the spectrum is consistent with PYTHIA |Sj01| (pQCD) calculations. 
Both direct reconstruction and indirect semi-leptonic decay electron measurements (spec- 
trum and elliptic flow) at top RHIC energy offer us an opportunity to understand not 
only pQCD, but also the early stage features of the matter created in heavy ion collisions. 



1.2.6 What have we learned so far? 



So what is happening after heavy ions collide? Fig. 11.131 is a cartoon showing the evo- 
lution of excitation energy E* (represented as temperature T for a thermally equilibrated 
system) and the particles being created and decoupling from the system. There is no 
absolute scale on the plot because we don't know most of these variables yet. 
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Figure 1.13: A time scale plot to show the evolution of the created matter. The several 
ellipses depict the system shape during evolution. And creation and decoupling of 
different identified particles are also shown on the plot. 



At RHIC, after two heavy nuclei collide, plenty of energy (~ 70% of the initial en- 
ergy per participant |Bea04j ) is deposited into a compact overlap region suddenly. The 
initial Bjorken energy density is well above the critical energy density for a phase tran- 
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sition |Adl04aj . Those particles from hard processes (such as jets and heavy quarks) as 
well as direct photons {qg —>■ q'-f, qq g'~f) are created in these initial collisions. The 
overlap region is an ellipse in coordinate space. Then rescatterrings between components 
in the system make the system expand and approach thermal equilibrium. Soft mini-jets 
{Pt <~ 2 GeV/c) are created abundantly. The anisotropy in coordinate space quenches 
due to the expansion, but the anisotropy in momentum space saturates quickly during 
this stage |KH03j . The system will sooner or later reach its maximum temperature, and 
the matter at this stage is expected to have the partonic d.o.f. equation of state. Fast 
partons traversing through this hot dense matter will lose large amounts of their energy 
(jet quenching). Meanwhile, particles (photons, leptons, g, u, d, s, but very small frac- 
tion of heavy quarks) with relative small masses are produced thermally in this stage. 
With the continuous expansion of the system, the temperature starts to drop, the matter 
becomes more and more dilute and those partonic components begin to hadronize from 
the system either through coalescence or fragmentation. In the stage after hadronization, 
hadrons keep rescatterring between each other and fractions between different hadrons are 
still varying. During this stage, some hadrons with small scattering cross sections start 
to freeze out from the system with fixed momentum such as charm hadrons {D, J/ip), 
multi-strange hadrons {Q, 0, S). Then the system becomes dilute enough that there is no 
inelastic scattering between different hadrons and the fractions of different hadrons are 
stable. This point is called chemical freeze-out. Some hadrons (vr, K, p{p)) in the system 
keep scattering elastically and finally all particles freeze out from the system with each 
particle having stable momentum. This is called kinetic freeze-out. 

In heavy ion collisions, what we want to study ultimately is the equilibrated matter 
with partonic d.o.f.. Multi-strange hadrons, and more early freeze-out charm hadrons 
are ideal probes to illustrate the properties of this matter created in collisions. Based on 
the knowledge we have learned from three years of RHIC runs, we need these probes to 
answer the questions: whether the partonic matter is locally thermalized or not, and if 
yes, what are its thermal and symmetry properties? 
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CHAPTER 2 



Experimental Set-up 

2.1 RHIC accelerator 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
(BNL) is designed to accelerate and collide heavy ions and polarized protons with high 
luminosity, allowing physicists to explore the strong interaction through many extensive 
and intensive measurements. It is the first facility to collide heavy ion beams and the top 
center-of-mass collision energy is 200 GeV per nucleon pair, which is about more than 10 
times greater than the highest energy reached at previous fixed target experiments. The 
purpose of this extraordinary new accelerator is to seek out and explore new high-energy 
forms of matter and thus continue the centuries-old quest to understand the nature and 
origins of matter at its most basic level. RHIC is also delivering polarized proton beams 
up to center-of-mass energy 500 GeV to carry on vigorous spin scientific programs. 

Fig. 12.11 shows a diagram of the RHIC machine complex, including a Van de Graaff 
facility, a linear proton accelerator, the booster synchrotron, the Alternative Gradient 
Synchrotron (AGS), and ultimately the RHIC synchrotron ring. For Au beam operations, 
the Au ions with charge Q = —le are created using the Pulsed Sputter Ion Source. Then 
they are accelerated through the Van de Graaff facility and a series of stripping foils, and 
the Au ions at the exit are with a kinetic energy of 1 MeV/nucleon and a net charge 
of Q = +32e. The ions are then injected into the booster synchrotron and accelerated 
to an energy of 95 MeV/nucleon. After the Au ions leave the booster, they are further 
stripped to Q = +77e and transferred into the AGS, where they are accelerated to 8.86 
GeV/nucleon and sorted into four final bunches. Finally, the ions are injected into RHIC 
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and stripped to the bare charge state ofQ = +79e during the transfer. For p+p operations, 
protons are injected from the 200 MeV Linac into the booster, followed by acceleration 
in the AGS and injection into RHIC. 
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Figure 2.1: Schematic of the RHIC complex. RHIC's two 3.8-kilometer rings collide 
relativistic heavy ions and polarized protons at six intersection points. 



RHIC consists of two concentric super-conducting storage rings that are called blue 
and yellow rings, respectively. Each ring has its own dependent set of bending and focusing 
magnets as well as ratio frequency cavities, but both share a common horizontal plane 
in the tunnel. The rings have six interaction points, and 4 of them are equipped with 
detectors. They are two large experiments STAR (6 o'clock), PHENIX (8 o'clock) and 
two small ones PHOBOS (10 o'clock) and BRAHMS (2 o'clock), respectively. 

To date, RHIC has been run in p+p , d + Au , Au + Au and Cu + Cu configurations. 



20 



2.2 STAR detector 



The Solenoidal Tracker at RHIC (STAR) is a specially designed detector to track 
thousands of particles simultaneously produced by each ion collision at RHIC. It has 
an azimuthal symmetric acceptance and covers large range around mid-rapidity. STAR 
consists of several subsystems and a main tracker - the Time Projection Chamber (TPC) 
located in a homogenous solenoidal analyzing magnet. 



Magnet 




Figure 2.2: Cutaway view of the STAR detector. It includes the partial installed 
ElectroMagnetic Calorimeter (EMC) and two prototypes Time-of-Flight (TOP) 

detectors. 

Fig. l2.2l shows the cutaway view of the STAR detector. The main tracker - TPC covers 
the |?7| < 1.5 and 27r in azimuth. The details of TPC detector will be discussed in the next 
section. There are inner detectors Silicon Vertex Tracker (SVT) and Silicon Strip Detector 
(SSD) close to the beam pipe, which provides additional high precision space points on 
track so that it improves the position resolution and allows us to reconstruct the secondary 
vertex of weak decay particles. There are two Forward TPC (FTPC) detectors covering 
2.8 < \ri\ < 3.8 to track particles at forward and backward rapidity. One prototype 
tray of Time- Of- Flight (TOP) detector using scintillator materials (TOPp) was installed 
since Run II and another prototype tray of TOP detector using Multigap Resistive Plate 
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Chamber (MRPC) technology (TOFr) was installed since Run III. Each replaces one 
Central Trigger Barrel (CTB) tray (120 in total) surrounding TPC. They are used to 
test the performance of upcoming full barrel TOP detector upgrade which is expected to 
extend the PID capability of STAR greatly. Part of barrel ElectronMagnetic Calorimeter 
(EMC) was also installed since Run II. They are used to measure the electromagnetic 
probes - electrons and photons. 

There are some main trigger detectors: Zero Degree Calorimeter (ZDC), CTB and 
Beam-Beam Counters (BBC). Two ZDCs locates on each side ~ 18 m away from the 
collision points. Each is centered at 0° and covers ~ 2.5 mrad. The ZDCs are hadronic 
calorimeters to detect the outgoing neutrons. They are put beyond the dipole magnets 
which bend away the charged fragments. The ZDC signals are used for monitoring the 
heavy ion beam luminosity and for the experiments triggers. The CTB is a collection of 
scintillating tiles covering the whole barrel ektexine of the TPC. The CTB will be mostly 
used to select central triggered events in heavy ion collisions by measuring the occupancy 
of those CTB slats. The BBC subsystem covers 3.3 < \r]\ < 5.0, measuring the "beam- 
jets" at high rapidity from Non-Singly Diffractive (NSD) inelastic p + p interactions. It 
consists of two disk shaped scintillating detectors, with one placed at each endcap of the 
TPC (3.5 m from TPC center). Each BBC disk is composed of scintillating tiles that are 
arranged in a hexagonal closest packing. The p+p NSD trigger sums the output of all tiles 
on each BBC and requires a coincidence of both BBC's firing above noise threshold within 
a time window. Some other detectors are used for special triggers, e.g. pseudo Vertex 
Position Detectors (pVPDs) are used for TOP triggered events (this will be discussed in 
the following section) , and EMC is used to trigger on high pt particle events etc. 

The STAR magnet is cylindrical in design with a length of 6.85 m and has inner and 
outer diameters of 5.27 m and 7.32 m, respectively. It generates a field along the length 
of the cyhnder having a maximum of \Bz\ — 0.5 T. It allows the tracking detectors to 
measure the helical trajectory of charged particles to get their momenta. To date, the 
STAR magnet has been run in full field, reversed full field and half filed configurations. 
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2.3 Main tracker - TPC 



TPC is the "heart" of the STAR detector |Andn8j . Consisting of a 4.2 m long cyhnder 
with 4.0 m in diameter, it is the largest single TPC in the world. The cylinder is concentric 
with the beam pipe, and the inner and outer radii of the active volume are 0.5 m and 2.0 m, 
respectively. It can measure charged particles within momentum 0.15 < pt/ (GeV/c)< 30 
(0.075 GeV/c low limit for 0.25 T). The TPC covers the full region of azimuth (0 < < 
27r) and covers the pseodurapidity range of |?7| < 2 for inner radius and \ri\ < 1 for outer 
radius. Fig. 12.31 shows a cutaway view of the structure of the TPC. 




Figure 2.3: Cutaway view of the TPC detector at STAR. 



The TPC is divided into two parts by the central membrane. It is typically held at 
28 kV high voltage. A chain of 183 resistors and equipotential rings along the inner and 
outer field cage create a uniform drift filed (~ 135 V/cm) from the central membrane to 
the ground planes where anode wires and pad planes are organized into 12 sectors for 
each sub-volume of the TPC. The working gas of the TPC is two gas mixture — PIO (Ar 
90% + CH4 10%) regulated at 2 mbar above the atmospheric pressure. The electron drift 
velocity in PIO is relatively fast, ~ 5.45 cm/yus at 130 V/cm drift field. The gas mixture 
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must satisfy multiple requirements and the gas gains are ~ 3770 and ~ 1230 for the 
inner and outer sectors working at normal anode voltages (1170 V for inner and 1390 
V for outer), respectively. Each readout plane is instrumented with a thin Multi-Wire 
Proportional Chamber (MWPC) together with a pad chamber readout. Each pad plane 
is also divided into inner and outer sub-sectors, while the inner sub-sector is designed 
to handle high track density near collision vertex. 136,608 readout pads provide {x, y) 
coordinate information, while z coordinate is provided by 512 time buckets and the drift 
velocity. Typical resolution is ~ 0.5 — 1.0 mm. 

When charged particles traverse the TPC, they liberate the electrons from the TPC 
gas due to the ionization energy loss (dE/dx). These electrons are drifted towards the 
end cap planes of the TPC. There the signal induced on a readout pad is amphfied and 
integrated by a circuit containing a pre-amplifier and a shaper. Then it is digitalized and 
then transmitted over a set of optical fibers to STAR Data AcQuisition system (DAQ). 

The TPC reconstruction process begins by the 3D coordinate space points finding. 
This step results in a collection of points reported in global Cartesian coordinates. The 
Timing Projection chamber Tracker (TPT) algorithm is then used to reconstruct tracks 
by hehcal trajectory fit. The resulted track collection from the TPC is combined with 
any other available tracking detector reconstruction results and then refit by application 
of a Kalman filter routine — a complete and robust statistical treatment. The primary 
collision vertex is then reconstructed from these global tracks and a refit on these tracks 
with the distance of closest approach (dca) less the 3 cm is preformed by a constrained 
Kalman fit that forces the track to originate from the primary vertex. The primary 
vertex resolution is ~ 350 /im with more than 1000 tracks. The refit results are stored 
as primary tracks collection in the container. The reconstruction efficiency including the 
detector acceptance for primary tracks depends on the particle type, track quality cuts, 
Pt , track multiplicity etc. The typical value for the primary pions with Nf^ > 24 and 
\rj\ < 0.7, dca < 3.0 cm is approximate constant at px > 0.4 GeV/c: >~ 90% for Au + 
Au peripheral collisions and ~ 80% for central collisions, respectively. 

The TPC can also identify particles by the dE/dx of charged particles traversing the 
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TPC gas. The mean rate of dE/dx is given by the Bethe-Bloch equation ITTl |Eidn4j : 



dE 
dx 




max 



2 



6 



(2.1) 



The meaning of each symbol can be referred to |Eidn4j . Different types of particles 
(different rest masses) with the same momentum have different kinematic variables (3 
(7), which may result in distinguishable dE/dx. The typical resolution of dE/dx in Au 
+ Au collisions is ~ 8%, which makes the tt/K separation up to p ~ 0.7 GeV/c and 
proton/meson separation up to p ~ 1.1 GeV/c. 

A new recent technique was developed to identify high momentum (p > 3 GeV/c) 
pions and protons in the relativistic rising region of dE/dx |Xu04j benefiting from the 
advantage that the mean rates of dE/dx for different particles have a visible separation 
in the relativistic rising region (~ 2a separation for pions and protons). Due to large 
acceptance of the TPC, using the topology of their weak decay in the TPC, the K^, A(A) 
etc. can be identified across pt region 0.3 < px / (GeY /c)< 7.0 (upper edge limited by 
statistics). Resonances (-ft'*, 0, A etc.) can be reconstructed through the event mixing 
technique |Adl02bj . Fig. 12.41 shows the PID capabilities up to date with the TPC. In 
addition, the TOF PID capability is also shown on the plot which will be discussed in the 
next section and the analysis part as well. 

2.4 Prototype TOF detector 

STAR has proposed the full barrel Time-Of-Flight (TOF) detector upgrade based on 
the Multigap Resistive Plate Chamber (MRPC) technology in the coming future. The 
TOFp detector (a prototype based on scintillator technology) was installed since Run 
II |Llo04j . It replaced one of CTB trays, covering — 1 < < 0, and 7r/30 in azimuth. It 
contains 41 scintillator slats with the signal read out by Photo Multiplier Tubes (PMTs). 
The resolution of TOFp is ~ 85 ps in Au + Au collisions. However, due to the significant 
higher cost by the PMTs, this design will not be used in the full TOF upgrade. 
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Baryons 
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Figure 2.4: reach of particle identification capability with STAR detectors for Run 
II and Run III configurations. The upper edges of rdE/dx, weak decay topology, event 

mixing are limited by statistics. 

In Run III and Run IV, new prototypes of TOF detector based on MRPC (TOFr) 
were installed. Each also replaced one CTB tray, covering — 1 < < and 7r/30 in 
azimuth too. In Run III, 28 MRPC modules were installed in the tray and 12 of them 
were equipped with electronics, corresponding to ~ 0.3% of the TPC acceptance [AdaOScj . 
In Run IV, 24 modules were installed in a new tray and the tray was put in the same 
position in STAR as Run III (but slightly global z position shift), but only 12 modules 
were equipped with valid electronics, which means the acceptance in Run IV was roughly 
similar to that in Run III. 

Two pVPDs were installed as well since Run II to provide a starting time for TOF 
detectors, each staying 5.4 m away from the TPC center along the beam line |Llo04j . 
Each pVPD consists of three detecting element tubes covering ~ 19% of the total solid 
angle in 4.43 < \ri\ < 4.94. Due to different multiplicities, the effective timing resolution 
of total starting time is 25 ps, 85 ps and 140 ps for 200 GeV Au + Au , d + Au and p + p 
collisions, respectively. 

MRPC technology was first developed by the CERN ALICE group. Fig. 12.51 shows 
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Figure 2.5: Two-side view of a MRPC mo dule |Shan2j . 



Table 2.1: TOF system performance in Run III and Run IV 



Timing Resolution (ps) 





pVPDs 


TOFr system 


TOFp system 


Run 


+ Au @ 200 GeV 


85 


120 (85)^ 


100-140 


III 


p + p @ 200 GeV 


140 


160 (85) 


N/A 




Au 


+ Au @ 62.4 GeV 


55 


105 (89) 


110 (95) 


Run 


Au + Au 


FF/RFF, w/o E pVPD^ 


40 


95 (86) 


96 (87) 


IV 


@ 


FF/RFF 


27 


86 (82) 


92 (88) 




200 GeV 


HF-^ 


20 


82 (80) 


85 (83) 



^ numbers in the brackets are those for intrinsic TOFr/TOFp detector. 
^ east pVPD is not available for those runs from 5022001 to 5036000. 
* HF - half field; FF - full field; RFF - reversed full field. 
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0.5 1 1.5 2 2.5 3 3.5 

Particle momentum p (GeV/c) 

Figure 2.6: 1/13 vs. momentum (p) from 200 GeV d + Au collisions. Separations 
between pions and kaons, protons and mesons are achieved up to px ~ 1-6 and ~ 3.0 
GeV/c, respectively. The insert plot shows th? = p'^{l/j3^ — 1) for 1.2 < < 1.4 Gev/c. 

the two side views (long edge view on top and short edge view on bottom) of an MRPC 
module appropriate for STAR |Sha02j . An MRPC basically consists a stack of resistive 
plates with a series of uniform gas gaps. It works in avalanch mode. Electrodes are 
applied to the outer surface of the outer plates. With a strong electric field applied on, 
the internal plates are left electrically floating and they will keep the correct voltage due 
to the flow of electrons and ions created in avalanches. There are six read-out strips on 
each module in this design. The first beam test for 6-gap MRPCs at CERN PS-TIO 
facility with piab = 7 GeV/c pions beam resulted in a ~ 65 ps timing resolution with more 
than 95% detecting efficiency and the module is capable of working at high event rate 
(500 Hz/cm2) |Sha02j . These modules were then assembled in a prototype TOP tray and 
tested in the AGS radiation area. Similar resolution was obtained. In RHIC Run III and 
Run IV, the MRPC modules in TOFr trays installed in the STAR detector were applied 
on the high voltage of 14 kV and with the working gas of 95% freon and 5% isobutane. 
The charged particle detecting efficiency is > 95% at high voltage plateau. 

TOP system calibrations include the start time calibration from pVPDs and TOPr/TOPp 
flight time calibration. The main sources need to be considered are global time offset due 
to different electronics delays, the correlation between the amplitude and the timing 
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signals, the correlation between the hit position and the timing signals etc. Detailed 
calibrations on TOF systems can be found in |Rua04bl IAda03cj (TOFr) and |Llo04j 
(TOFp). Table. 12.11 lists the calibrated timing resolution results for TOFr and TOFp 
system in Run III and Run IV (Run IV 200 GeV Au + Au results are from test sample 
so far). The results show that the intrinsic timing resolution of TOFr was ~ 85 ps and 
this performance was quite stable in two-year runs. Fig. 12.61 shows the hadron PID capa- 
bility of TOFr system in 200 GeV d + Au collisions |Ada03cj . The performance of TOFr 
satisfied STAR TOF system upgrade requirements. 
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CHAPTER 3 



Single electron transverse momentum distributions 

Single lepton measurement has been proposed as a feasible and effective way to extract 
heavy flavor production since long time ago. Since the decay kinematics of heavy flavor 
hadrons to leptons is well known, the transverse momentum (px) distribution of single 
lepton can reveal that of heavy flavor hadrons. And furthermore, many topics may be 
carried on, such as total cross section, heavy quark fragmentation, heavy quark energy 
loss in medium in A + A collisions, etc. In this chapter, analysis details of single electron 
Pt distributions from 200 GeV d + Au , p + p collisions and 62.4 GeV Au + Au collisions 
will be presented. 

3.1 Single electrons from d + Au and p + p collisions at s/snn = 
200 GeV 

3.1.1 Data sets and Trigger 

In RHIC Run III, STAR detector has accumulated several data sets from d + Au 
and p + p coUisions at ^/sN^ = 200 GeV. Besides the minimum bias triggered data from 
d + Au collisions and the Non Singly Diffractive (NSD) collisions from p + p collisions, 
special triggers were set up for TOFr detector data accumulation in both d + Au and 
p + p collisions since TOFr detector has a relative very small acceptance. Table ITT] lists 
all the data sets under the selections used in this analysis from Run III. 

The d + Au minimum bias trigger required at least one beam-rapidity neutron in 
the ZDC in the Au beam outgoing direction, which is assigned negative pseudorapid- 
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Table 3.1: Data sets from Run III used in this analysis 



Collision System 


Trigger 


Vertex Z Selection 


Data Sample Size 


d + An 


minbias 




Vz 


< 30 cm 


6.42 M 




dAuTOF 




Vz 


< 50 cm 


1.89 M 


p + p 


NSD 




Vz 


< 30 cm 


4.35 M 




ppFPDTOFu 




Vz 


< 50 cm 


1.08 M 



ity (r/) |Ada03bj . The trigger accepts 95 ± 3% of the d + Au hadronic cross section. 
The NSD p + p events are triggered on the coincidence of two BBCs, which are angu- 
lar scintillator detectors covering 3.3 < |?7| < 5.0 |Ada03dj . The NSD cross section was 
measured to 30.0 ± 3.5 mb by a van der Meer scan and PYTHIA simulation of the BBC 
acceptance |Ada03cj . The TOF trigger set up in Run III was to select events with a valid 



pVPD coincidence and at least one TOFr hit (out of 72). The trigger enhancement factors 
in comparison to minimum bias trigger are ~10 in d + An and ~40 in p + p collisions. 

Centrality definition ind + An collisions was based on the charged particle multiplicity 
in —3.8 < 1] < —2.8, measured by the FTPC in the Au beam direction |Ada03bj . But in 
this analysis, due to the limited statistics, the analysis was done on minimum bias d + 
Au events, not on each specified centralities. 

3.1.2 Electron identification and hadron contamination 

TPC is the main detector in STAR for tracking and identifying charged particles. 
With TPC only, electron identification is complicated because the electron band crosses 
the hadrons bands. A previous analysis on conversion electron was done for Run II Au + 
Au data using pure topological method, however, electron identification was limited up 
to pt ~ 0.8 GeV/c jAda04d[ rfoh02| . 

In addition to its hadron identification capabilities, TOFr detector can be used to 
identify electrons in combination with the dE/dx in the TPC, shown in Fig. 13.11 The 
top panel of left plot shows the 2-D scattering plot of dE/dx vs. particle momentum {p) 
for the charged particles with good TOFr matched hits from d + Au collisions |Rua04bj . 
The selection criteria are listed in Table 13.21 If additional particle velocity {(3) cut 
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— 1| < 0.03 is applied on, this plot is shown in the bottom panel. Slow hadrons 
were eliminated and electrons band, then, can be separated from those fast hadrons band 
more readily. With combination of dE/dx in the sTPC and j3 in TOFr, electrons can be 
identified starting from p ~ 0.15 GeV, while the upper limit in momentum is constrained 
by the statistics in this analysis. 




0.5 1 1.5 2 2.5 



Particle momentum p (GeV/c) 

Figure 3.1: Left: Electron identification by combining (3 from TOFr and dE/dx from 
TPC With the cut of — 1| < 0.03, electron band can be separated from the 
hadrons readily, shown in the bottom panel. Right: The dE/dx projection plot in 

1.0 <pT/(GeV/c)< 1.5. 

Electrons were selected under the following cut: 

dE/dx{p) > 2.4 + 0.65 x (1 - e-(P-o-i^)/o-i) + 0.1 x p (3.1) 

where p is in the unit of GeV and dE/dx is in keV/cm. Hadron contamination under 
this cut was studied by fitting to the dE/dx distributions in each bin. The 2-Gaussian 
function cannot describe the shoulder region of fast hadron peak very well in lower prp 
region, a function of exponential+gaussian was used in the fit. At ~ 2 — 3 GeV, 
statistics cannot enable us to distinguish the difference of these two fitting methods, so 
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2-Gaussian fit was performed in tliis pt region. Fig. 13. 21 shows the fitting results in several 
typical Pt bins from d + Au collisions. The arrows denote the cut from Eq. 13.11 Hadron 
contamination ratio was estimated from these fits, shown in Fig. 13. 4t and the efficiency 
for electrons was corrected for in the final spectra and will be discussed in the sections 
later. 
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Figure 3.2: dE/dx projections in several pt bins with the combination of TOFr PID 
cut. The electron yields were extracted from the fit. The green arrows denote the 
dE/dx cut to select electrons in each pt bin. 



At higher pt (2 — 4 GeV/c), electrons could be identified directly in the TPC since 
hadrons have lower dE/dx due to the relativistic rise of the dE/dx for electrons. Positrons 
are more difficult to be identified using dE / dx alone because of the large background from 
the deuteron band. Additional track quality cuts were used in this selection: number of 
fit points is required to be at least 25 and the number of dE/dx points is required to be at 
least 16. Fig. 13. 31 shows the logarithm of dE/dx distributions in each pt bin. A 3-gaussian 
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function (e + vr + K/p) fit to the distribution was performed based on the assumptions that 
all the particles have the same dE/dx resolution and one gaussian function can describe 
kaons and protons. The results are shown in Fig. 13.31 Electrons with dE/dx greater than 
the peak mean value were selected. 
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o 


0.03413 ± 0.00015 


p_yield 


3.952e+04±315 


p_mean 


0.3882 ± 0.0003 


e_yield 


577.3 ± 58.6 


e_mean 


0.5701 ±0.0031 







Il7 / \ e/Ti(c;,>0)= 3.11 
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Figure 3.3: Logarithm of dE/dx distributions for negative charged particles in three pt 
bins between 2 — 4 GeV/c. A 3-Gaussian fit was used to extracted the electron yield and 
hadron contamination. The pink arrows denote the electrons selection cut. 

Hadron contamination was extracted from the fits in these two methods, respectively, 
shown in Fig. 13.41 for both d + Au and p + p collisions. The hadron contamination not 
only need to be subtracted statistically in the final yield, but also should be considered 
carefully on the track-by-track analysis later on. This will be discussed in detail in the 
next section. 

Low hadron contamination and high electron efficiency are two competiting aspects. 
We performed the cut scanning to optimize the electron selection. Fig. 13.51 shows the 
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Figure 3.4: Hadron contamination fractions as a function of Pt from two electron PID 
methods in ci + Au (left) and p + p (right) collisions. 

electron efficiency vs. hadron contamination. Since we are limited for the statistics 
currently, optimization results in the selection shown in Eq. 13.11 with ~ 90% electron 
efficiency and ~ 10% hadron contamination. Fig. 13.61 shows the pt dependent ratio of 
electrons to the total hadrons. It shows ~ 1% ratio at pr > I GeV/c. Combining the 
~ 10% contamination ratio under the selection shown in Fig. 13.41 the estimated hadron 
rejection power under such selection is ~ 10^ at pt > 1 GeV/c and ~ 10^ — 10^ at px ~ 0.5 
GeV/c. 
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Figure 3.5: Electron efficiency vs. 
hadron contamination fraction by 
varying the dE/dx cut to selection 
electrons for 4 p^ bins. 
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Figure 3.6: The ratio of electron yield 
to total hadron yield. By combining 
Fig. 13. 4t the hadron rejection power 
under the current selection (Eq. 13. 1|) 
was estimated to be ~ 10^ at > 1 
GeV/c. 
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Raw yields of inclusive electrons were then calculated from the fits to dE/dx distri- 
butions at each pt bin. Table. 13.21 lists the event-wise and track-wise cuts used in the 
electron selection. 



Table 3.2: Electron selection criteria 



Method 


TOF+rfE/rfx 


dE/dx 


triggerword 


2300 for d + ku 






1300 ioT p + p 




\VertexZ\ < 


50 cm 


30 cm 


primary track ? 


Yes 


Yes 


nFitPts ^ 


15 


25 


nFitPts/nMax > 


0.52 


N/A 


ndEdxPts ^ 


N/A 


16 


rapidity 


(-0.5, 0) 


(-0.5, 0) 


j3 from TOFr 


|l//3- 1| < 0.03 


N/A 


TOFr hit quality 


10 < ADC < 300 


N/A 




-2.7 < ziocai/cm < 3.4 


N/A 



3.1.3 Acceptance, Efficiency and Trigger bias 

Raw yields need to be corrected for the detector acceptance, reconstruction and se- 
lection efficiency etc., which were done through Monte Carlo (MC) simulations. Full 
detector simulations and/or embedding study show a ~ 90% reconstruction efficiency for 
e~^/e~ tracks, shown in Fig. 13.71 As for the tracks hitting TOFr detector, there are two 
ways to determine the acceptance/efficiency: the first one is divided into two steps, TPC 
reconstruction efficiency and the matching efficiency from TPC to TOFr from real data; 
the other one is get the total efficiency (including acceptance) from embedding data. De- 
tailed matching efficiency study is referred to |Rua04bj . Due to the lower statistics of 
electrons, matching efficiencies of vr^ got from real data was used for electrons and a cor- 
rection factor 90% was used to account for the loss (scattering and decay) of pions flying 
from TPC to TOFr. In the later method, d + Au HIJING events with an additional MC 
e~^/e~' track within TOFr acceptance were embedded into real zero-biased events, and the 
total efficiency was extracted from the embedding data directly. Fig. 13.81 shows the total 
efficiency (including acceptance) correction factor from two methods in d + Au collisions. 



36 



The efficiency correction in p + p collisions was used as the same number as those from d 
+ Au results because of the similar multiplicities. 
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Figure 3.7: TPC track reconstruction 
efficiency for electrons from embedding 
data. 
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Figure 3.8: Acceptance included 
efficiency for TOFr + dE/dx electrons 
selection from two methods. 



The additional TOFr PID cut — 1| < 0.03 used for slow hadrons rejection will 
naturally lead to some loss of the electrons. A simple fast simulation based on the TOFr 
timing resolution (120 ps in (i + Au and 160 ps in p + p ), charged particle phase space 
{Pt,v) distributions was used to determine this correction. The result shows ~ 95% 
efficiency in d + Au and ~ 87% efficiency in p + p under this cut and this correction was 
applied in addition to the others above. 

As for the TOF triggered data sets, the bias study of this trigger is mandatory. De- 
tailed investigation was introduced in |Rua04bj . The results show no significant pt 
dependence for the spectrum at px > 0.3 GeV/c according to the TOFr trigger and the 
charged particle multiplicity bias correction was done for both d + Au and p + p data 
sets. 

To improve statistics, the r] cut —0.5 < rj < was removed. This leads to a ~ 44% 
increase for the TOFr acceptance and no pt dependence was found. Then the yields will 
be corrected to a definite t] range for calculation according to this factor. The final results 
from loosing t] cut are consistent within statistical errors with those with rj cut applied. 
In the following text, we will remove this cut in the discussion. 
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3.1.4 Photonic background contribution 



There are several sources that can contribute to the final single electron spectra in 
STAR environment: 

• Photon (from tt*^, rj etc.) conversions in inner detectors 

• TT^, 7] etc. scalar meson Dalitz decay 

• p, uj, (f) vector meson di-electron decay and/or Dalitz decay 

• K"^ decay — Kes 

• heavy flavor (c, b) hadron semi-leptonic decay 

• other sources (Dell-Yan, heavy quarkonium decay, thermal electrons, direct photon 
conversion etc.) 

In this analysis, the first four sources are considered to be the background — photonic 
background, which need to be subtracted from the inclusive yield. As for the last item, 
theoretical predictions ||Sj01t [RapOH lAlaOlj show no significant contributions with large 
uncertainties. This part of contribution is neglected in this analysis. The heavy fiavor 
hadron decay is considered to be the only signal. 

From the previous measurement |Adc02j and the estimation from the STAR envi- 
ronment, electrons from photon conversion and vr'' Dalitz decay dominate the total yield, 
especially at low pt ■ At pt > 2 GeV/c, the signal may become visible due to the different 
shape between the total spectrum and background spectrum. In PHENIX collaboration 
single electron analysis, they used cocktail method to subtract each contribution one by 
Adc02j . However, the pr distribution and the total yield of each source particles 



one 



were used as assumed input in their analysis, which may introduce large uncertainties. In 
STAR, due to the large acceptance of TPC, background from photon conversion etc. may 
be reconstructed experimentally. The pioneering analysis using topological method were 
done in |Ada04dl IJoh02j . In the analysis we report here, since we have already tagged 
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one electron/positron track, we don't want to use the whole topological reconstruction. 
Instead, we try to use the kinematical features of photon conversion and 7r° Dalitz pro- 
cesses to justify whether the tagged electron/positron is from photonic background or not. 
Firstly, use TOFr or dE/dx at high pt to tag one electron/ positron track from the pri- 
mary collision vertex, then loop all other opposite-sign global tracks, whose helixes have 
a distance of closest approach (dca) to the tagged tracks' less than 3 cm, to find the other 
partner track in the TPC Because the TPC acceptance is large, the pair reconstruction 
efficiency is reasonable. 

To study the kinematics of the photon conversion and 7r° Dalitz decay processes, 
simulations of HIJING |WG91j for d + Au system and PYTHIA Sj01| for p + p system 
plus full detector description GEANT were investigated in detail. The data sets used 
are 1.38 M + Au HIJING (vl.35 + GCALOR) minimum bias events and 1.2 M p + p 
PYTHIA (v6.203 + GCALOR) events. The tagged electrons were required with the same 
cut as those for data, but all good TPC tracks were selected to increase the statistics. 
Fig. 13.91 shows the invariant mass square and opening angle distribution of electron pairs 
from different sources. 
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Figure 3.9: Invariant mass square and opening angle distributions of electron pairs 
from several sources in simulation. The blue arrows depict the cut thresholds in this 

analysis. 



As expected, the results show those background mostly has small invariant mass square 
and/or small opening angle. This offers us a good opportunity to estimate the single elec- 
tron background with almost topology-blind method. Those tagged electrons/positrons 
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with the invariant mass square or opening angle below certain thresholds will be consid- 
ered to be background from photon conversion and/or 7r° Dalitz decay, and be rejected 
in the analysis. 

To understand the background reconstruction capability in our environment, simula- 
tion study were done on the same data sets. Fig. I3.1(JI shows the reconstruction efficiency 
of photon conversion and 7r° Dalitz decay in the TPC The simulation shows an almost 
constant efficiency ~ 60% at pt above 1 GeV/c. From the vr^ embedding data, in which 
MC 7r° tracks through full detector response simulation were embedded into real zero-bias 
events , the similar result within errors was also obtained. 
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Figure 3.10: Photon conversion and 
Dalitz decay background 
reconstruction efficiency from 
simulations. 



Figure 3.11: The ratio of single 
electrons from photon conversion to 
those from vr*^ Dalitz decay from 
simulations. 



By applying this efficiency to the spectrum after the kinematic cuts, one can extract 
the estimated total background from the main sources: photon conversion and 7r° Dalitz 
decay. Other contributions as well as the fractions of those could not be extracted from 
data directly, and we determined them from simulations. Fig. 13.111 shows the ratio of 
electrons from photon conversion to those from 7r° Dalitz decay. The electron spectrum 
from vr*^ Dalitz was then used as a reference for those contributions from other light meson 
decays. Fig. 13.121 shows the spectra of electrons from different light hadron decays and 
the contribution of each source w.r.t 7r° contribution from d + Au HIJING simulations. 
From this , we found that the total contribution of photon conversion and 7r° Dalitz 
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decay constitutes ~ 95%, which were measured experimentally. Table 13.31 shows the 
fractions of each electron background sources to the total background, averaging over 
1 < ]9'jr'/(GeV/c) < 3. And the upper panel of Fig. 13.131 shows the inclusive electron 
spectra and the photonic background contributions in d + Au and p + p collisions. Clear 
excess of electrons at > 1 GeV/c is shown and the bottom panels of that plot shows 
the ratio of inclusive electrons to the total background electrons. The grey bands depict 
a 20% systematic uncertainty of the total background. 
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Figure 3.12: The single electron spectra from various source contributions from 
simulations. The bottom panel shows the relative fractions to 7r° Dalitz decay, which 
were used in the estimation of other source contributions in data analysis. 



3.1.5 Signal extraction 

One simple way to extract the signal so far is to just subtract the estimated total 
background from the total inclusive electron spectrum. In practice, to reduce the statis- 
tically subtracting errors, we subtracted the background from the left spectrum after the 
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Pt (GeV/c) p-r (GeV/c) 

Figure 3.13: Upper panels: Electron distributions from p+p (left) and d+Au (right) 
collisions. Solid and open symbols depict electrons/positrons ((e^+e^)/2) identified via 
a combination of TOFr and dE/dx, and electrons (e~) identified via dE/dx alone. The 
measured total photonic backgrounds are shown as solid lines. Dashed lines depict the 
various contributing sources. The fractions were derived from simulations. Bottom 
panels: the ratio of inclusive electrons to the total backgrounds. The gray band 
represents the systematic uncertainty in each panel. 



Table 3.3: Photonic sources contributions to single electron spectrum 



Source 


fraction (%) 


7 conversion 


80 ± 7 


7r° Dalitz 


15 ± 2 


V 


3.2 


CO 


0.99 




0.22 


P 


0.19 


Ke3 


0.21 
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kinematical cuts applied. Besides the real background one can reject under such cuts, 
some signals may also be cut out due to the random combination with the hadron tracks 
(mostly pions), which need to be compensated after subtracting those background from 
photon conversion and 7r° Dalitz decay. Hadrons contaminating in the electron selection 
will be considered as following the random rejection power behavior under the above kine- 
matical cuts. Fig. 13.141 shows the random rejection fractions as a function of pt from the 
MC studies. 
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Figure 3.14: The electron fraction rejected by random particles as a function of pt 

The signal fraction was then extracted from the following equations: 



eff X B + ran x {S + H) = cut 
;i - eff) X 5 + (1 - ran) x {S + H) = sur 

{S + B)/H = r 



(3.2a) 
(3.2b) 
(3.2c) 



where B, S, H represent background electrons, signal electrons and contaminating 
hadrons, respectively, eff denotes the efficiency of the background reconstruction, and 
ran denotes the random rejection efficiency for signals and hadrons. cut and sur represent 
the yields being cut out and surviving cuts in data, respectively. The electrons to hadrons 
ratio (related to hadron contaminations) r has already been obtained from previous fit to 
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dE/dx distributions. From ()3.2|) . we have 

(e// + ran/r) x sur — (1 — e// + (1 — ran)/r) x cut 

{eff - ran) x (1 + 1/r) ^ 

when hadron contamination is neghgible ( r ^ 1 ), the above formula can be written as 

g ^ e// X sur - (1 - e//) x cwt ^ 
e// — ran 

Signal yields were extracted in each pt bin. 

The bin-by-bin systematic errors are dominated by the uncertainties of the background 
reconstruction efficiency from MC simulations and the hadron contamination. The errors 
were then calculated by varying the efficiency to Icr uncertainty away and changing the 
dE/dx cut in the electron selection for different hadron contamination levels. Table EiU 
list the non-photonic electron spectra data points and the errors. 

Table 3.4: Non-photonic electron spectra, the numbers represent yields ± statistical 

errors ± systematical errors. 







PT (GeV/c) 


data points 






1.0 


- 1.2 


(1.30 ±0.71 ±0.83) X 10- 


'A 






1.2 


- 1.4 


(5.52 ±3.66 ±2.41) x 10" 


-5 




TOFr + 


1.4 


- 1.6 


(6.42 ±2.43 ±3.40) x 10" 


-5 


d + Kn 


dE/dx 


1.6 


- 1.8 


(3.88 ±1.61 ±0.88) X 10" 


-5 


minbias 




1.8 


-2.0 


(1.55 ±0.85 ±0.26) x 10" 


-5 






2.0 


-3.0 


(4.82 ± 1.99 ±0.47) x 10" 


-6 






2.0 


-2.4 


(7.53 ± 1.12 ±2.41) X 10" 


-6 




dE/dx 


2.4 


-3.0 


(1.19 ±0.41 ±0.36) X 10- 


-6 






3.0 


-4.0 


(2.45 ±1.22 ±1.72) x 10^ 


-7 






1.0 


- 1.2 


(2.23 ± 1.36 ±1.43) x 10^ 


~5 




TOFr + 


1.2 


- 1.6 


(7.02 ±4.51 ±3.30) x 10" 


-6 




dE/dx 


1.6 


-2.0 


(3.59 ±2.05 ±0.79) x 10" 


-6 


p + p 




2.0 


-3.0 


(7.16 ±4.57 ±0.70) x 10" 


'7 


NSD 




2.0 


-2.4 


(1.10 ±0.30 ±0.36) X 10- 


-ti 




dE/dx 


2.4 


-3.0 


(1.89 ±0.85 ±0.68) x 10" 


-7 






3.0 


-4.0 


(4.11 ±3.02 ±1.68) X 10- 


-8 
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3.1.6 Background subtraction check 



The photon conversions from the inner detectors are dominated by the inner SVT 
detector and its supporting materials. These materials are mainly within \z\ <~ 50 cm 
|,Toh02j ■ As a double check the background subtraction, we made the vertex z position 
distributions of electron tracks normalized to the event vertex z position distributions, 
shown as solid symbols in Fig. l3.15l upper panel. The "M" structure reveals the conversion 
material positions. The open symbols in the same panel shows the same distributions for 
the cut out electrons, which are considered as conversion background. Both distributions 
show similar structure. If we take into account the background efficiency estimated from 
the discussion above, and subtract them from the total distribution, the result is depicted 
in the bottom panel in that figure. We cannot see significant material effect. The flat 
distribution demonstrates that the background subtracting method is reasonable. 
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Figure 3.15: Normalized vertex z position distributions of electron tracks for total, 

cut-out background and expected signal. 



Another important check is the electron/positron symmetry. Fig. 13.161 shows the non- 
photonic yield of electrons and positron separately. Within statistical errors, they are 
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consistent with each other. This indicates the possible 6 electrons , which are knocked 
out from atoms in material by charged particles, does not contribute significantly in the 
final spectrum. 




0.5 1 1.5 2 2.5 

Transverse momentum (GeV/c) 



Figure 3.16: Non-photonic electrons and positrons spectra from d + Au collisions. 



3.1.7 Electron dN/dy and implications for total charm production cross sec- 
tion 

Many previous measurements tried to extract the charm production cross section from 
electron spectra only |Bot90| IAdc02j . However, this method is strongly model dependent 
for the charm hadron spectrum is unknown and as well as the fractions between 
different charm hadrons. In this analysis, the electron spectrum only covers 1—4 GeV/c, 
which is corresponding to ~ 2 — 6 GeV/c in the charm hadron spectrum. The direct 
extraction from electron spectrum will have large systematic uncertainties. 

On the other hand, STAR collaboration has measured charm mesons through hadronic 
decay channels m. d + Au collisions jRua04a| ITai04j . The low pt D^{D^) (we use 
implying [D^ + D^)/2 in the following text) from channel D^{D'^)—>- K^ix^ [Branching 
ratio (B.R.) = 3.8%) was reconstructed through event mixing method |Adin2bj . In the 



following part, we will perform a fit combining both and electron data points in d + 
Au collisions. 
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Firstly, the decay kinematics of charm hadrons to electrons should be fixed. Two 
different decay packages were used: one is the GENBOD from CERNLIB |CERj . which 
does the multi-body decay in phase space; the other one is the PYIENT function from 
PYTHIA |Sj01| , which is believed to be a more reliable one since other effects (such as 
spin etc.) have been taken into account. We implement both for the systematic studies. 
Fig. 13.171 shows a 2D plot of the correlation between of charm hadrons and decayed 
daughter electron px- The result shows the parent charm hadrons have much higher p^ 
in comparison with the decayed electrons. 




1 2 3 4 5 

Pt (D) (GeV/c) 



Figure 3.17: 2D correlation plot of pt for charm hadrons and decayed electrons. The 
{Pt) of parent charm hadrons for decayed electrons at pr = 2 — 3 GeV/c is ~ 4.8 GeV/c. 



The fractions of different charm hadrons are also needed for the decay. Previous 
measurements as well as the model predictions show almost consistent ratios in broad 
energies and different collisions systems |Tai04j . We quote the values from e~^e~ collisions 
at = 91 GeV from PDG |Eidn4j : R = Ndo/N^c = 0.54 ± 0.05. And ^ 0.4 



from STAR d + Au measurement. According to the large uncertainty of the B.R. of 
—>■ + X and A+ e++X, and both are similar to that of D^, and the distributions 
of decay daughter electrons have similar shapes, we took and A^ as Thus the 
final input charm mesons are: 79% and 21% D+. This ratio will be tested for the 
systematic study later. 
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Table 3.5: The Combined fit characters for and electrons. 



par. 


scan range 


step 


result 


xVndf 


dN/dy{D^) 


(0.02, 0.04) 


0.0004 


0.029 ±0.004 




(Pt) 


(0.9, 1.6) 


0.05 


1.15 


18.8/9 


n 


(7, 16) 


0.5 


11 





The combined fit was then performed under the following assumptions: 

• charm meson spectrum follows a power law function up to ~ 6 GeV/c. 

• similar px spectrum shape between different charm hadrons. 

• the signal electrons are all from charm decays (bottom quark contribution will be 
discussed in chapter 5). 



The power law function is written as the following: 



_ 2(n- l)(n-2) dN 

2'npTdpTdy 7r(n - 3)2(pj.)^ dy 



Pt 



n — 3 



\ 



(Pt) 



(3.5) 



There are only 3 free parameters: dN/dy, (px) and the power n. A specified charm hadron 
with a set of these parameters was input into the generators, and the decayed electrons 
spectrum was obtained. A 3-dimensional scan on the {dN/dy, (pr), n) "plane" was done 
to fit and electron data points simultaneously. The point with the smallest value 
was set to be the fit result, was calculated from the following equation. 



X 



D 



O-D 



(3.6) 



where yn, ye denote the measured yields of and electrons and ctd, de denote the 
measured errors, /d, /e denote the expected values from input power law function for 
and electrons respectively. To avoid the px position issue in large px bins, we used the 
yield dN instead of dN/pxdpT in each px bin. Table EiH lists all the fitting characters. 

The error estimation was through the contour scan in the 3-D "plane" with the = 
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X,r, 



+ 1. The error of dN/dy was then obtained by projecting this 3-D contour into 



dN/dy axis. Fig. 13.181 shows the plots of corrected non-photonic electrons spectra and 
the combining fit results for and electrons spectra in c? + Au collisions. The curve for 
p + p collisions is by scaling down that of + Au collisions by {Nun) = 7.5 ± 0.4. 
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Figure 3.18: Reconstructed (solid squares) px distributions from d+Au collisions 
and non-photonic electron pt distributions from p+p collisions (triangles) and d+Au 
collisions (circles). Solid and dashed lines are the fit results from both and electron 
spectra in d+Au collisions. The dotted line is scaled down by a factor of 
{Nbin) = 7.5 ± 0.4 from d+Au to p+p collisions. 



Once the dN/ dy was extracted, the total charm production cross section per nucleon- 
nucleon interaction at RHIC energy can be calculated from Eq. 13.71 



cr. 



NN 



da 



NN 



dy 



y=0 



xf 



dNoo 
dy 



=0 X i? X 



cr; 



pp 

inel 



{Nbin) 



(3.7) 



In this equation, the factor R is the fraction in total charm hadrons, as mentioned 
before, and / is a factor when extrapolating the dN/dy at mid-rapidity to full rapidity. 
This factor / was extracted from PYTHIA model study and we estimated a ~ 15% 
systematic error on this factor from different parameters in that model. The charm 
production cross section per nucleon-nucleon collision at mid-rapidity and total cross 
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section at y/s = 200 GeV is: 



da 



NN 



dy 



{y/s = 200 GeV) = 0.30 ± 0.04 mb 



(3.8) 



a^,^{Vs = 200 GeV) = 1.42 ± 0.20 mb (3.9) 

3.1.8 Systematic error study 

The systematic errors contributing to the final total charm cross section were studied. 
One of the important sources is the uncertainty of charm hadron pt spectrum and the 
fractions of different charm hadrons. The effect caused by the charm hadron pt spectrum 
was studied by changing the power law parameters to la away. The correlations between 
parameters were not taken into account in the estimation, instead, we looped all the 
combinations to find the largest deviation as the estimation systematic uncertainty. The 
effect caused by the uncertain fractions of different charm hadrons was studied by changing 
the effective fraction in the decay : D+ = 0.79 : 0.21 to 0.90 : 0.10 and 0.70 : 0.30. 
For the charm decaying to electron kinematics, the difference of two decay codes was take 
into the systematic errors. Table 13.61 shows the contributions from possible sources. 

Table 3.6: Systematic errors to the final charm cross sections. 



Source 


Relative contribution 


charm hadron pt spectrum 




13% 


fractions of charm hadrons 




10% 


decay codes 




12% 


NFitHits 15 ^ 25 




14% 


rapidity distribution 




15% 


{Nbin) in ci + Au 




5% 


Normalization 




10% 


Sum 


31% 



50 



3.2 Single electrons from Au + Au collisions at s/snn = 62.4 
GeV 

The purpose of the data analysis oi d + Au and p+p colhsions is to set up the basehne 
for heavy ion colhsions. Since the large data sample of Au + Au collisions at ^/SNN = 
200 GeV is not available yet, the smaller sample of Au + Au collisions at ^/SNN = 62.4 
GeV has been analyzed to gain understandings of the necessary techniques for Au + Au 
collisions. Due to the significant increase of multiplicities in Au + Au collisions, and the 
consequent decrease of track qualities, the situation in Au + Au collisions becomes more 
complicated. In this section, I would discuss the inclusive electron spectrum from TOF + 
dE/dx method and the photonic background estimation. Since the charm yield is pretty 
low compared to photonic contributions, we would not expect to extract the signal with 
acceptable errors. 

3.2.1 Data set and electron PID 

RHIC has offered a relative short time Au + Au beams at ^/snn = 62.4 GeV and 
STAR has accumulated ~ 15 M events in total. With the minimum bias trigger (0 — 80%) 
and vertex z position selection, the useful physics events number in this analysis is ~ 6.4 
M. 

The calibration results for TOF detectors are ~ 110 ps resolutions for both TOFr and 
TOFp system, with ~ 55 ps start timing resolution included. The hadron PID capability 
was reported in jSha05j . Similar as before, electrons can be identified by combining TOF 
and dE/dx in the TPC. Primary electrons were selected under the criteria shown in 
Table. IT71 

Fig. 13. 191 shows the dE/dx vs. particle momentum after a (3 cut from TOF. Electrons 
band can be separated from hadrons. The dE/dx resolution in Au + Au collisions is 
worse than that in d + Au and p + p due to much higher multiplicities, so the separation 
in this plot is not as good as that in d + Au and p + p , see Fig. 13.11 We performed both 
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Table 3.7: Electron selection criteria in Au + Au 



Method 


TOF+dE/dx 




Vertex Z\ < 


30 cm 


primary track ? 


Yes 




nFitPts ^ 


25 


ndEdxPts ^ 


15 




rapidity 


(-1.0, 0) 




X^/ndf 


(0., 3.0) 


p from TOF 


\l/f3-l\ < 0.03 


TOFr hit quality 


30 < ADC < 300 






-2.7 < ziocai/cm< 3.4 






\yiocai - yd < 1-9 cm 


TOFp hit quality 


thi < ADC < th2 






2.0 < ziocai/cm< 18.0 






0.4 < yiocai/cm.< 3.2 



two-gaussian function and exponential+gaussian function fit to extract the electron raw 
yields in each pt bin. 
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Figure 3.19: dE/dx vs. particle momentum after a TOF (3 cut {\l/(3 — 1| < 0.03) 

Hadron contamination becomes larger if we select electrons with the same efficiency 
as before. We need to optimize the pure electron sample selection for the background 
study later. Fig. 13.201 shows the hadron contamination fractions under different electron 
selections. 

The fit results show that to have contamination under control (< 10%), we need the 
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Figure 3.20: Hadron contamination fractions for different electron selections in Au + 
Au 62.4 GeV. Statistics are poor at > 2 GeV/c. 

selection (jg > with electron efficiency ~ 50%. 

3.2.2 Non-photonic background estimation in Au + Au 

In Au + Au collision events, due to large multiplicity, the PID-blinding partner track 
reconstruction will lead to a huge combinatorial background. Even a critical dE/dx 
selection is applied, this combinatorial background is still significant. Fig. 13.211 shows 
the electron pair candidate invariant mass distribution. The tagged electron track was 
selected from TOF with the cuts shown in Table. 13.71 and additional < cTg < 3 (we will 
open it to —1 < cTg < 3 later). The partner track candidate was selected under the criteria 
shown in Table. 13.81 

Table 3.8: Partner candidate selection criteria in Au + Au 



charge 


opposite to tagged track 


primary/global ? 


global 


nFitPts ^ 


15 


nFitPts/nMax > 

X Indf 


0.52 


(0., 3.0) 




(-1., 3.0)t 


dca of e+,e^ 


(0.0, 3.0) cm 



several different cuts were tried. 



TOF Au+Au @ 62.4 GeV 
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Figure 3.21: Invariant mass of electron candidates pair from the same event in Au + 
Au 62.4 GeV collisions. The plot shows significant combinatorial contribution near the 

conversion peak. 

To reconstruct the photon conversion and tt*^ Dalitz, we need to understand the com- 
binatory background shape in Fig. 13.211 In this step, the photon conversion and other 
correlations will be considered as signals, and combinatorial background was produced 
by rotating the partner track momentum 'p' — "p', then do the normalization of the 
mass spectrum in the region 0.8 < Me+e-/(GeV/c^)< 2.0. Fig. 13.221 shows this results in 
several px bins. Here pr is the px of tagged electrons. The plots show the combinatorial 
background was very well reproduced. In the first plot, there are two small bumps at 
Me+e- ~ 0.3 GeV/c^ and 0.5 GeV/c^. These are due to the misidentified pions which 
come from Kg and p, u decays, respectively. This has been confirmed through MC decay 
simulations. 

The photon conversion peak is clearly seen near zero mass region, and the offset from 
is due to the opening angle resolution |Joh02j . The 7r° Dalitz contribution is not visible 
in this case, possibly because the Dalitz contribution is much smaller and distribution 
is much broader compared to conversion processes. We subtracted the combinatorial 
background from the real distribution, and integrated the remaining distribution from 
— 0.15 GeV/c^ to get the reconstructed main photonic background raw yield. In this 
case, we assumed that both photon conversion and vr'' Dalitz decays were reconstructed. 
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Figure 3.22: Main pliotonic background reconstruction in Au + Au collisions. The 
crosses depict the real electron pair candidate invariant mass distributions and the 
histograms represent the combinatorial backgrounds. The photon conversion peak is 
clear in ~ mass region, while the 7r° Dalitz is hard to see. 
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The background raw yield need to be corrected for the reconstruction efficiency, as we 
did it for d + Au and p + p colhsions. This efficiency was calculated from Au + Au 62.4 
GeV HIJING events plus full detector MC simulations. After \Vz\ < 30 cm cut, ~ 53 K 
events were used in the calculation. We took all TPC electron tracks without any dE/dx 
and TOF hit cut to improve the statistics. The procedure is the same as we did before. 
Fig. 13.231 shows the background reconstruction efficiency in Au + Au 62.4 GeV compared 
with d + Au results. Due to the higher multiplicities, the relative lower reconstruction 
efficiency does make sense. 

— \ ' 1 ' 1 ' \ ' \ ' 

O d+Au 200 GeV HIJING 

• Au+Au 62.4 GeV HIJING 

:©:::g:::^::i::^^:::::::i::::::::4::::: 

I \ I I I I I \ I \ I 

0.5 1 1.5 2 2.5 3 

Transverse momentum pj (GeV/c) 

Figure 3.23: Photonic background reconstruction efficiency from Au + Au 62.4 GeV 
HIJING simulations. Also shown on the plot is that from d + Au 200 GeV HIJING 

simulations. 

This efficiency was used to correct for the photonic background raw yield obtained 
from above. In addition, a ~ 5% fraction of other photonic background (from d + Au 
results) and the cXe selection efficiency were also included. Then we can get the corrected 
background spectrum compared with the inclusive spectrum, shown as Fig. 13.241 All these 
were not corrected for the single track efficiency and acceptance yet. The reconstructed 
background matches the total inclusive electron spectrum. This is not surprised since the 
charm yield is pretty low at 62.4 GeV. 

This improved method for Au + Au system working makes us confident to extract the 
charm signal at 2 — 3 GeV/c from the coming 30 M minimum bias 200 GeV Au + Au 
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Figure 3.24: Raw inclusive electron yield and the reconstruction efficiency corrected 
photonic background contributions under different electron/positron track selections. 
The numbers in the brackets on the plot show the cTg cut for the tagged track and the 



3.2.3 Spectra from Au + Au collisions at ^/snn = 62.4 GeV 

Single electron track efficiency and acceptance is needed to correct the raw yield spec- 
trum. The efficiency calculation is similar to what has been introduced before. Here, 
the method combining TPC tracking efficiency from MC simulation and TOF matching 
efficiency from real data was used. Fig. 13.251 shows the matching efficiency of pions from 
TPC to TOF (including both TOFr and TOFp) within — 1 < y < window |Sha05j . 

The TPC tracking efficiency from d + Au was lowered down by 5% to as an estimation 
for Au + Au 62.4 GeV collisions according to the multiplicity difference in d + Au and 
Au + Au collisions. The TOF velocity selection efficiency was estimated to be ~ 95%. 
With all these efficiency corrected, we can get the spectrum of inclusive electrons, shown 
in Fig. 13.261 Also shown in the plot are inclusive charged pion spectrum from TOFr 
measurement in Au + Au 62.4 GeV and a previous ISR non-photonic electron spectrum 
from p + p 52.7 GeV collisions scaled by the Nbin to Au + Au collisions. This plot shows 
the expected non-photonic signal is about an order of magnitude lower than the inclusive 



partner track, respectively. 



data. 
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Figure 3.25: Matching efficiency (inciuding detector response) from TPC to TOF 
(including botfi TOFr and TOFp) of pions wfiin — 1 < y < window. 

electron spectrum. This is also consistent with the photonic background reconstruction 
discussed in the previous section. 

Although more than ~ 90% of the electrons up to 3 GeV/c are photonic background 
in Au + Au 62.4 GeV so that we cannot extract the non-photonic signal with reasonable 
errors, this photonic background reconstruction method works well and technically, it is 
ready for the coming 200 GeV large data sample analysis to answer some of the issues 
about the charm production in heavy ion collisions. 



58 



10^ 



> 

O 



"D , 
Q. 



10' 



▲ 

o 



7t* @62GeV 

inclusive (e*+e')/2 @62GeV 
non-photonic (e*+e')/2 @52.7GeV- 
p+p scaled to AuAu 



O 



O 



o 



o 



o 



Transverse momentum (GeV/c) 



Figure 3.26: Inclusive electron spectrum in Au + Au 62.4 GcV compared with the 
charged pion spectrum in the same system and the ISR non-photonic electron spectrum 
in 52.7 GeV p + p collisions scaled with A^^m to Au + Au minimum bias collisions. 
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CHAPTER 4 



Single electron azimuthal anisotropy distributions 

The data set used in this analysis is Au + Au minimum bias triggered (0 — 80%) 
events at ^/snn = 62.4 GeV. This data sample is the same as that described before. 
The anisotropic parameters were obtained using the event plane analysis technique. The 
detailed analysis method introduction can be found in |PV98j . 



4.1 Event plane and its resolution 

In heavy ion collisions, the event plane is reconstructed from the detected final particle 
azimuths. The acceptance and efficiency of the detectors in azimuth was corrected by 
compensating the azimuth to a fiat distribution with (j) weights. Technically, the (p weights 
were created for different days to deal with the different situations in a long running 
period. Additional pt weights were also applied to improve the event plane resolution. 
The second order harmonic event plane azimuth can be calculated from the Q vector, 
as Eq. OlOl 

^2 = ( arctan ^ j /2, < ^2 < vr (4.1) 

Q = {Q^,Qy) = (^Wi-cos{2(l)i), ^Wi-sin(2(/.i)j (4.2) 

Here, Wi is the weight for each track included in the event plane calculation, which includes 
both the (f) weight and the pt weight. Fig. 14. II shows a set of typical weight distributions 
for a single day. The bumps at ~ 6 in "East" and "Fast East" plots are due to a bad 
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sector on the east side of the TPC The tracks selected in the event plane calculation 
should satisfy the criteria listed in Table. 14.11 
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Figure 4.1: weight used for a single day for minimum bias 62.4 GeV Au + Au 
collisions. The bump around ~ 6 rad in (f) in east and far east set is due to a bad sector 

in east side of TPC. 



Table 4.1: Track selection in event plane calculation 



nFitPts 


> 15 


nFitPts/nMax 


> 0.52 


Pt 


(0.1, 4.0) 


V 


(-1.3, 1.3) 


global dca 


(0.0, 3.0) 



Fig. 14.21 shows the second order event plane azimuthal angle distribution. This distri- 
bution was fit to a constant value and the fit quality /ndf = 126/119 indicates a good 
event plane reconstruction. 

The resolution of the event plane was calculated using the sub-event method |P V98j . 
Each event was divided into two random sub-events with nearly equal multiplicity. The 
event plane was reconstructed in each sub-event, denoted as \l/2 ^2- Then the event 
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plane resolution r = (cos [2(\l/2 — ^rp)]) can be calculated from Eq.(14) and (11) from 
[PV98| : 



(C0S[2(^2 - ^rp)]) 



X2exp(-x^/4)x[/o(x^/4) + /i(x^/4)] 



(4.3) 



(cos[2(vI/^ - vl/,^)]) = v'(cos[2(M/^ - ^)]) 



(4.4) 



X2 = ^2/0- = V2V2N 



(4.5) 
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Figure 4.2: The second harmonic event plane azimuthal angle distribution. A constant 
fit with /ndf close to unity means a good event plane reconstruction. 

Firstly, we obtained the sub-event resolution (cos[2(\l/2 — ^rp)]) from Eq. 14.41 Then 
Eq. 14.31 can be solved as an iterative routine to extract the sub-event X2- This variable is 
proportional to \fN according to Eq. l4.5| so the total event X2 was obtained by X2 — ^^X2- 
After putting this X2 i^^^ 14. 3( we calculated the final full event resolution. The 
physical V2 is calculated as V2 = v^'^ /r, where v^'^ is the observed V2 and r is the event 
plane resolution. However, experimentally, what we observe is (^2^''), the averaged v'^'^ 
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over a data sample, such as the minimum bias events. Then 



,obs 



{but ^) 



(r) 



(4.6) 



This is not quite correct when we just divide the {v2^^) by the event-wise averaged 
resolution (r). Practically, we calculated a track-wise averaged resolution by weighting 
each event with the number of observed particles in the V2^^ calculation |Sor04aj . Fig. 14.31 
shows the event-wise averaged and track-wise averaged event plane resolution for 9 cen- 
trality bins and minimum bias (0 — 80%) events. The final event plane resolution used 
for correction is 68% for minimum bias Au + Au 62.4 GeV collisions. 
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Figure 4.3: Event plane resolutions for all centrality bins (bin 1-9) and minimum bias 
(0 — 80%) events (bin 0). Black points show the resolution with event weight only. Red 
points show the improved resolution with multiplicity weight. 



4.2 Elliptic flow of stable hadrons 

The TOF detectors were used for hadron identification in this analysis. A multiple 
gaussian fit was applied to the distributions in each px bin. In those pt bins where 
peaks start to merge, a cut of more than 95% purity was applied in the particle selec- 
tion. In the pt region beyond the tt and K separation with TOF only, dE/dx was used 
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additionally to help identify pions while K PID is limited by its low yield |Sha05j . 

There is a certain correlation between the selected particle (f)j and the event plane 
azimuth calculated including this particle. This so-called auto correlation was removed 
by excluding this selected particle in a new event plane azimuth calculation. The new Qj 
vector was constructed by: 



Qj = {Qjx, Qjy) = i^^i- cos(20j), ^ Wi ■ sin(20. 



(4.7) 



-1 Qjy 



^2, = (^tan-1 ^ J /2 (4.8) 
Then the f 2'''' in each pj- bin was calculated by 

vt = {cos[2{<jy,-^2,)]) (4.9) 

Fig. 14. 41 shows the measured resolution corrected V2 of vr^, and p{p) from the TOF 
detector for Au + Au 62.4 GeV minimum bias (0 — 80%) collisions. The results show: 
at low pt region, the mass ordering is similar to what has been observed in 200 GeV Au 
+ Au collisions; and at intermediate pt { > 2 GeV/c ) the proton t>2 overtakes pion t>2, 
which is also similar to that in 200 GeV case. To check the NCQ scaling that was observed 
in 200 GeV data, a function parametrization Eq. 14.101 was used to fit to the measured 

1 — exp[—[pT/n — b)/c\ 
where n is the number of constituent quarks. 

The fit results are shown in Table. 14.21 The consistency of the fit parameters for 
and p + p indicates an agreement with NCQ scaling. This scaling in 62.4 GeV was also 
further demonstrated by the measurements of V2 of Kg, A, S and fl jSor04bj . 

Comparisons of the identified particle f 2 from different energy measurements are shown 
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Figure 4.4: Elliptic flow V2 of vr^, and p{p) from TOF detector at Au + Au 62.4 
GeV minimum bias (0 — 80%) collisions. Curves depict the parametrized function fit 

from |DESn4j . 
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Figure 4.5: Minimum bias Au + Au 62.4 GeV V2 (data points) compared with those 
from 200 GeV measurements in Au + Au collisions at RHIC (non-solid lines) j Adl03aj 
and 17.2 GeV measurements for pions and protons (solid lines) in Pb + Pb collisions at 

SPS [Alt03j . 
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Table 4.2: Parametrization for V2 of and p + p 







p + p 


a 


6.3 ± 1.8 


6.1 ±1.3 


b 


0.44 ±0.04 


0.38 ±0.03 


c 


0.14 ±0.06 


0.13 ±0.04 


d 


0.05 ± 1.05 


0.66 ±0.92 



in Fig. 14.51 Those are measurements from ^JsNN = 200 GeV Au ± Au minimum bias 
(0 — 92%) collisions jAdl03aj and ^/SNN = 17.2 GeV Pb ± Pb minimum bias colli- 
sions jAlt03j . The results show the V2 from 62.4 GeV are very similar to those from 
200 GeV measurement, while they are significantly different from those in 17.2 GeV mea- 
surement. 

Charged pions can also be identified using the relativistic dE/dx {idE/dx) method |-X-u04j 
up to ~ 7 GeV/c (limited by the statistics). A simple selection of < 0-7^ < 2 was used to 
select a pion sample and the contamination from kaons and protons was estimated to be 
less than 3% pCHoi] . Fig. lOl shows the pion f 2 measured over a large px region ( 0.2 — 7.0 
GeV/c ). The event plane method shows a continuous fiat region up to almost 7 GeV/c. 
However, the non flow effect in the high px region becomes important and even more 
significant in central and peripheral collisions jAda04ej . To answer the question of high 
Pt identified particle f2, we need other methods (such as multiparticle cumulant methods 
etc jAda04f] ) to minimize the non flow effects. 



4.3 Elliptic flow of inclusive and photonic background electrons 

The electron PID capability in Au ± Au 62.4 GeV was already shown in Fig. 13.191 
With the selection of < de < 3, we identify electrons with the purity more than 90% in 
1 — 3 GeV/c. The inclusive electron V2 was extracted using a method similar to that used 
for V2 of vr^, and p{p). Fig. 14.71 shows the results compared with V2 of vr^. Since the 
dominant source of the electrons are photon conversion where the photons are from two 
photon vr'^ decay processes, and vr'' Dalitz decays, the electrons at lower px will carry the 
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Figure 4.6: Pion V2 from TOF, TOF + dEdx and relativistic dE/dx {idE/dx) 

measurements in Au + Au 62.4 GeV. 

anisotropy of those parent pions at higher p-p. This will be discussed in the next chapter. 
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Figure 4.7: V2 of inclusive electrons and pions from TOF detector in Au + Au 62.4 

GeV minimum bias (0 — 80%) collisions. 



Photonic electron V2 is needed to extract the interesting non-photonic electron V2- 
We used the event mixing technique to reconstruct the main photonic background from 
the TPC only. The tagged and partner electron candidates are selected using the same 
selection criteria shown in Table. 13.71 and 13.81 but we used all TPC candidates instead 
of TOF tracks only. Both tagged and partner track are required to have < cTg < 3 to 
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reject hadrons. The event buffer used for mixing is divided into a 10 x 10 x 9 lattice in 
{V:,, Mult, ^2) with lattice size of 6 cm x 60 x 20°. Two 3-D {M^+e~, pr , ^<P = (p - ^2) 
histograms for the same event and mixed-event were stored. The additional dimension 
in \1'2 for event mixing is to control the residual background in different bins. The 
histograms were then projected onto the M^+e- axis for each {px, A0) bin. The mixed 
event distribution, which is expected to describe the combinatorial contribution, was 
normalized to the same event distribution in the high M^+e- region, where there are 
no significant decay correlations. The normalized mixed event distribution was then 
subtracted from the same event distribution. The photonic electrons were extracted for 
every A0 bin in each pt bin. The A0 distribution was fit to the following function to get 
the observed fn''*. 



ocl + 2t;f*cos(2A0) (4.11) 



d{A(f)) 

Fig. 14.81 shows the Mg+e- distribution of the electron candidate pairs from the same 
event (crosses) and the mixed event (histograms) in different A0 bin for tagged electrons 
2.0 < pr/(GeV/c)< 2.5. The combinatorial background shows oscillation in A0, which 
is expected from the anisotropy from final state hadrons. After mixed combinatorial 
background subtracted, the photonic source electron anisotropy is shown in Fig. 14.91 Since 
the photon conversion (tt^ Dalitz continuum may hide in) distribution is not a gaussian- 
like peak, it is hard to do a signal-|-background fit as we did for other resonances to extract 
the signal yields. We need to understand the residual background well. So in the yield 
extraction in each A0 bin, we tried several residual background estimations to extract the 
systematic errors: (i) no residual, (ii) use a constant fit in 0.15 < Me+e-/(GeV/c^)< 0.25 
and take it as a residual, (iii) the same as (ii), but take the averaged value through all A0 
bins as a common residual. Fig. 14.101 shows the yields from each A0 bin using method 
(iii). Eg. 14. Ill was used to fit the distribution and the result is ^2''* = (13.4 ±4.4)%. After 
resolution correction, we can get the V2 of electrons from photonic sources. Fig. 14. Ill shows 
the photonic electron f 2 as a function of pt, compared with the inclusive electron f 2 and 
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pion V2- The photonic electron V2 is consistent with that of inclusive electrons, which is 
expected from the spectra analysis in the previous chapter. 

In general, the V2 of inclusive electrons is the sum of photonic (B) and non-photonic 
(S) electron t>2, weighted by each yield fraction fs, fs- 

vT = fB^ + fsx (4.12) 
So the , the elliptic flow of non-photonic electrons, can be calculated as 

^ _2 JB 2_ ^4_^3^ 

Js 

From above formula, the ratio of S/B is needed in this analysis. So the spectra analysis 
is necessary for the V2 of non-photonic electrons. 

A recent technique was developed to deal with such signal f 2 extraction from back- 
ground |Sor04aj . We can get the V2 of inclusive electrons in different Mg+e- bins — V2{M). 
And if we know the signal to background ration r(M), then we can use two free parameters 
ff and V2 to fit the V2{M). This statistical fit method looks promising in the V2 calcula- 
tion for Vo particles, and vr^, p,p in dE/dx relativistical rising region. So this method is 
feasible for the extraction of photonic electrons f 2 out of the mixture with combinatorial 
background, and also that of non-photonic electrons f 2 out of inclusive electrons. 
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Figure 4.8: Invariant mass Mg+e- distributions of electron pair candidates and the 
combinatorial background from rotating method in each A0 bin for tagged electron 

2.0 <pT/(GeV/c)< 2.5. 
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Figure 4.9: Invariant mass M^+e- distribution of combinatorial background subtracted 

electron pair candidates in each A0. 
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Figure 4.10: A0 distribution of photonic electrons in 2.0 < pT/{GeV/c)< 2.5. The 
distribution was fit to [0] x (1 + 2f2^'^' cos(2A(/))) to extract the observed f2'"''. 
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Figure 4.11: Inclusive electrons V2 (from TOP trays) and photonic electrons V2 (from 
the TPC) compared with charged vr^ V2. The grey bands in each pt bin for photonic 
electrons f 2 depict the systematic uncertainties from the combinatorial background 

estimation. 
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CHAPTER 5 



Discussion 

5.1 Open charm production in high energy coUisions 

5.1.1 Total charm cross section 

As we discussed in the introduction chapter, the heavy quark total cross section mea- 
surement offers a powerful test for pQCD calculations. Since plenty of measurements were 
made at low energies while few were made at ^/s > 100 GeV, theoretical calculations are 
often tuned to match the low energy data points and then extrapolated to high energies. 
Because the parameters in the calculations (scales, heavy quark mass etc.) are not under- 
stood well yet in all energies, the predictions at high energies differ significantly. Fig. 15.11 
shows almost all the total charm cross section measurements made so far compared with 
several theoretical calculations. 

The low energy {y/s <100 GeV) data points were taken from the review paper |Tav87j 
(which summarized all the charm measurements before 1987) and from several experi- 
ments after 1988 ( refer to |Adc02j for data points)^. A factor of 1.5 was applied on all 
these data points to include the contribution of Df and A^ |FMN97j . We do not include 
total cross sections from those references that were extrapolated from high xp and/or 
had extremely low efficiency from correlation measurements, (with extrapolation factor of 
^ 10). The PHENIX and UA2 experiments extracted the cross section from the single 
electron spectrum in 130 GeV Au + Au and 630 GeV p + p collisions, respectively. 

The lines depict different theoretical calculations and model fits. PYTHIA, a model 
^Detailed data points selections are listed in Appendix C. 
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Figure 5.1: Total cc cross section per nucleon-nucleon collision vs. the collision energy 

(a/s ). The dashed line depicts a PYTHIA calculation with a specific set of 
parameters SjOl . The solid and dot-dashed lines depict two NLO pQCD calculations 

with MRST HO, rric = 1.2 GeV/c^, fip = and specified shown on the 
plot |Vog02| . The dotted line (in pink) depicts a power law fit to the data points, with 

power n ~ 2. 

based on pQCD calculations, is often used for predictions. The dashed line depicts 
a typical pQCD calculation with the only change of PDF from the default PYTHIA 
(CTEQ5L^CTEQ5M1). The result underpredicts the measurements at RHIC. Although 
the default PYTHIA takes initial radiation into account, high order processes seem to 
still not be well predicted. These include initial and final radiation, gluon splitting, and 
production through parton showers etc.. Recent analysis for CDF energy open beauty 
measurements showed these processes contribute a large fraction in the heavy fiavor cre- 
ation at Tevatron |Bus04j . This discrepancy indicates these processes may also play an 
important role at RHIC energy. 

The two curves labeled "NLO pQCD" depict the full NLO pQCD calculations from |VogU2 
Vog04| . 'Full" means for both PDFs and differential cross sections. The red dot-dashed 
line shows a calculation tuned to the low energy data points. This calculation underpre- 
dicts our measurements at 200 GeV as well. The scale parameters used are rric = 1.2 
GeV/c^, fiF = fJ'R = 2mc. When the NLO pQCD calculation is carried out with fin = rric 
instead of 2mc, it seems to reproduce the measurement at 200 GeV, but misses other data 
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points. An even smaller scale (/ir = mc/2) was tried, but it overshoots our measurement 
by almost an order of magnitude. The K-factors in pQCD calculations also differ sig- 
nificantly when changing the scales |Vog021 |Vog04| . As we know, when the scales go to 
~ 600 MeV, which is close to Aqcd, the feasibihty of pQCD is doubtful. The discrepancy 
indicates these scales can be energy dependent. 

On the other hand, recent tt" spectrum measurements at both mid-rapidity |Adl03bj 
and forward rapidity |Ada04g| at RHIC seem to offer clear evidence that pQCD can 
reproduce the vr*^ spectrum very well. This agreement even reaches as low as ~ 1 
GeV/c. Because of the large Q^, one might naturally think the agreement for heavy flavor 
production should be better. However, our measurement indicates a negative result, which 
means there can be some differences in the calculations between light flavor hadrons and 
heavy flavor hadrons. 

We fit a power law function to the data points in Fig. 15. II and the result is depicted as 
a pink dotted line with the power ~ 1.9 ± 0.2. This power law dependence has already 
been proposed in |RRL03j and they gave the power n ~ 1.6. The power law dependence 
of charged hadron production [n ~ 0.3) has been observed in A + A collisions from AGS 
to RHIC ABR03J and also predicted in the saturation approach |KL01j . Also for pion 
production, a linear scaling as a function of F = {^fs^N — 2mMY^'^ / \J ^]vn ~ ^/y/^NN was 
observed |Afa02j in central A + A collisions, which is equivalent to a power law dependence 
of ^/sNN with n ~ 0.5. For the kaons that have strange valence quarks, similar scaling 
seems to hold from AGS to RHIC, but with a larger power value n ~ 0.6 — 1. Hence, this 
power scaling for charm quark production up to RHIC is not surprising. One interesting 
thing is that heavier quark production has a larger value of power, which may be caused 
by a threshold effect in the production: once the production channels are opened, or above 
and near production threshold, the total cross section distribution on different flavors is 
likely to favor heavier quarks. 

Fig. 15.21 shows da/dy compared to theoretical predictions with smaller systematic 
uncertainties. This figure shows clearly that the theoretical calculations underpredict out 
measurement. 
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Figure 5.2: da /dy of charm quarks from STAR measurement compared with different 
theoretical predictions |VogU4| ITTFoHl lijC^HOI-jj . 

5.1.2 Spectrum comparison 

Open charm hadrons, including D'^{K-tt+), D*+{D°n+), D+{K~n+n+) and D%K-n+p), 
were reconstructed directly in c? + Au collisions (the measurements were statistically lim- 
ited in p + p ). A consistency check between the open charm spectrum and its expected 
semi-leptonic decay electron spectrum were made. The PYTHIA function for particle 
decay was used for generating the electron spectrum. The input open charm spectrum 
was fit to a power law function, and the parameters are |Tai04j : dN/dy = 0.0265 ±0.0035, 
(Pt) = 1-32 ± 0.08, and n = 8.3 ± 1.2. Fig. 15.31 shows the electron spectra compared with 
the expected contribution from charm quark semi-leptonic decays. The bands include 
the uncertainties from those parameters. The comparison demonstrates a consistency 
between the two measurements. 

There are already several electron spectra from charm semileptonic decay measure- 
ments. In Fig. 15.41 they are shown. 

The comparison shows a continuous increase in the yields with the increase of beam 
energies. However, from RHIC 130 GeV to UA2 630 GeV, the errors on the data points 
conceal the detailed structure. Not only the total yield, but also the shape of the spectra 
seem to become harder with the increase of collision energy. 
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Figure 5.3: Non-photonic electron spectra compared with the expected ones from the 
semi-leptonic decays of measurement open charm spectrum. 




Figure 5.4: Non-photonic electron spectra from several collision energies and collision 
systems. Data points are taken from |B76t IAdc02t lAdlOSl IBotQOj . Curves are expected 
spectra from models or measured open charm spectrum. 
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Figure 5.5: Measured non-photonic electron spectrum compared witli PYTHIA 



Fig. 15.51 shows a comparison of our measurement to PYTHIA calculations using dif- 
ferent parameters. The PYTHIA calculation from the PHENIX 130 GeV Au + Au 
paper |Adc02j was used as a starting line (black solid curve). This calculation deviates 
significantly from the result oi d + Au measurement at ^/snn = 200 GeV: it shows a 
lower yield and a steeper pt spectrum. Several parameters in PYTHIA were tuned to try 
to match our measurement: we increased the abundance of parton showers, we included 
higher order processes and gluon splitting processes, we changed the PDFs, we increased 
the initial (kj-) and we changed the charm quark mass. Although one of the calculations 
(red solid line) seems to match the spectrum, the large initial {kj^) requirement is ques- 
tionable. Theoretical predictions |Vit03j as well as measurements show a much smaller 
(kx) is needed for the 200 GeV d + Au system. In general, PYTHIA results seem to be 
hard to match to our measured spectrum. 

In the above comparisons, we didn't change anything about the charm quark frag- 
mentation function, which is a Peterson function in the default PYTHIA |Sj01| . The 
bare charm quark spectrum from PYTHIA seems to match the open charm measure- 
ment, which means a very hard fragmentation function (almost a 6 function) is needed 
for charm quarks |Tai04j at RHIC energy. This may indicate that charm quark coales- 



calculations with different parameter settings at 200 GeV. 
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cence processes become competitive with fragmentation processes in the 200 GeV d + Au 
system. A recent study shows this final state effect may play an important role through 
soft and shower recombination in d + Au collisions |HY04| IRua04bj . 

5.1.3 Bottom contribution 

In the above discussions, the bottom contribution to the non-photonic electron spec- 
trum was neglected. Since there is no bottom measurement around RHIC energy, this 
can only be studied through models. We used the default PYTHIA as a first estimation. 
Fig. l5.6l shows the results from this estimation where we used the scale a^c '■ cr^^i = 400 : 1 as 
the normalization |Vog02[rkP03j . At pt above 1 GeV/c, all the electrons from bottom are 
from direct semi-leptonic decays. This contribution is ~ 25 — 40% for 2 < pj'/{GeV/c)< 3 
and ~ 60 — 100% for 3 < pT/{GeV/c)< 4 to the charm contribution. The corre- 
sponding bottom contribution to the total non-photonic electrons is ~ 20 — 30% for 
2 < pT/(GeV/c)< 3 and ~ 40 - 50% for 3 < pT/(GeV/c)< 4. If we subtract this esti- 
mated bottom contribution from our spectra, the change in the total cross section from a 
combined fit is ~ 0.09 mb, within the current statistical error of our result. This fraction 
is strongly model dependent however, so we didn't include the bottom contribution in the 
final calculation. 

We also tried several of the parameter settings from Fig. 15.61 to estimate the bottom 
contribution. Almost all these results show the bottom contribution will begin to overcome 
the charm contribution around ~ 3 — 4 GeV/c. However, in these calculations, the 
heavy quark fragmentation functions have not been modified. As we discussed above, 
there is a hint that heavy quark may have a hard fragmentation function. In this case the 
difference of the fragmentation between charm quarks and bottom quarks may be much 
smaller than the default Peterson settings in PYTHIA. A recent study shows that under 
such assumptions, the bottom contribution on the electron spectrum does not overcome 
the charm contribution until px ~ 10 GeV/c |Lin04j . Precise measurements of heavy 
flavor spectra are needed to address these open issues. 
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Figure 5.6: Comparison between tlie contributions from charm quarks and bottom 

quarks from the defauh PYTHIA. 

5.1.4 Cronin effect of charm hadrons in + Au 

The well-known " Cronin effect" was named after the observation that particle (differ- 
ential) cross sections depend on the atomic number (A) and the particle species from low 
energy p + A collisions jCro73j . It is also characterized by the spectrum broadening in 
the intermediate px in p + A collisions compared to p + p collisions. The initial projectile 
partons' multi-scattering with target partons has been proposed as an explanation for this 
effect jWanOOt IVit03j . However, a recent measurement on the pseudorapidity asymmetry 
of charged hadron spectra shows a contrary result with these model predictions. Another 
approach through the final state recombination of soft and shower partons seems to re- 
produce the identified particle spectra in d + Au collisions |Rua04b| IHY04j . Measuring 
the Cronin effect of charm hadrons offers an additional opportunity to study the initial 
{kx) broadening and coalescence of the final state charm quarks. Even though we don't 
have enough statistics to measure the charm hadron spectra in both d + Au and p + p 
collisions, estimates of R^au^ nuclear modification factor of non-photonic electrons can 
give us some hints about these processes for charm hadrons. 

Fig. l5.7l shows the R'^j^^ for 1 < pT/{GeV/c)< 4 compared with a theoretical prediction 
curve based on the initial multi parton scattering picture. We used the measured charm 
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Figure 5.7: Cronin effect of non-photonic electrons in d + Au collisions, compared with 

the theoretical prediction of charm hadrons. 

hadron spectrum in d + Au collisions and generated that in p + p collisions according to 
the RdAu fi'om |Vog03| . Then we obtained the decayed electron spectra from both and 
calculated the predicted R^au- -Due to poor statistics, the current measurement cannot 
give definite conclusions. The averaged R'^au in 1 < pT/(GeV/c)< 4 is 1.3 ± 0.3 ± 0.3, 
which is consistent with N^in scaling within present uncertainties. 



5.2 Closed charm production 

5.2.1 J/ip production in p + p collisions 

The heavy quarkonia production mechanism in hadron hadron collisions is not well 
understood. Understanding this production mechanism in p + p collisions is needed to 
search for J/ ip suppression as a signature of QGP formation in central Au + Au collisions 
at RHIC energy. Recent theoretical progress |NLC03j and the measurement from PHENIX 
collaboration jAdl04bj strongly support the Color Octet Model (COM) for J/tp production 



in p + p collisions at ^/s = 200 GeV. Fig. 15.81 shows the measured J/ ip production cross 
sections compared with the total charm quark pair production cross sections in different 
energies. The ratio of (Jji^/acc is on the order of 10~^ and shows a slight decrease from 
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low energies up to RHIC energy. 
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Figure 5.8: ajj^ as a function of collision energy in comparison to the total charm 
cross section. J production mechanism is proved to validate the COM model at 

RHIC energy. 

5.2.2 J/ip production in heavy ion collisions 

In the hot dense QGP matter, the color screening effect is predicted to prevent heavy 
quark pairs from forming quarkonium bound states, and thus leads to a decrease of the 
ratio of closed charm (bottom) to open charm (bottom) jMS86j . This heavy quarko- 
nium suppression has been proposed as a signature of QGP formation. In low energy 
SPS 17.2 GeV Pb + Pb collisions, the NA50 experiment reported a suppression of 
heavy quarkonium suppression relative to normal nuclear absorption. This was inter- 
preted as an anomalous suppression due to the dissociation of cc pairs interacting with 
co-movers |Abr01j . 

In high energy heavy ion collisions, since there may be multiple cc pairs produced 
in a single collision, charmonium production through coalescence of charm quarks may 
become important: D + D J/ip + X. Since cc pairs are mostly produced from initial 
scatterings jLin96j . the partition of these charm quarks between open and closed ones is 



82 



deduced in statistical coalescence models |ABR03j from Eq. 15.11 . 



Nl°' = ]^N,, + N,, (5.1a) 

Noc oc gcfioc (5.1b) 
Nccocglricc (5.1c) 



In these models, the total charm yield is needed as an input to calculate Qc- The models 
often use 300 — 400 /zb for at RHIC energy, which is a typical value from NLO pQCD 
calculations j MNR93t |Vog02[ lRP03j . With this input, the statistical coalescence model, 
without any additional absorption effects, seems to agree with the J ftp measurement in 
Au + Au collisions within errors jAdl04cj . However, if our measurement o"^^ = 1.4 mb is 
used as the input, assuming the statistical coalescence models work in central Au + Au 
collisions at RHIC and ignoring the absorption effect, one would expect at least a factor 
of 3 increase above the number of binary scaling in the J/ip yield in central Au + Au 
collisions. The upper limit from PHENIX measurement on J/ ip production in central Au 
+ Au collisions seems to invalidate this expectation. 

Additional improvements for the nucleon absorption effects on J/ip were taken into 
account in |GR01j . But in central collisions, the production through coalescence can still 



make up a significant portion of the total production. Fig. 15.91 shows the predictions from 
different calculations. 

From Fig. 15. 9[ because of the dominant contribution from charm quark recombination, 
the expected J/ip yield in central Au + Au collisions inferred from our d + Au result is 
still significantly higher than expectations from number of binary scaling. Even when 
nuclear effects in the d + Au system are taken into account, which means acc ~ 1.0 mb 
is used as the input, this prediction is still above number of binary scaling. In the QGP, 
due to color screening, the J/ip yield will be suppressed comparing to Nbm scaling. Since 
most of the suppression will occur at low pt , the (pt) of J/ ip will be efficiently increased. 
On the other hand, if the QGP is formed, partons will reach kinetic equilibrium |ABR03j 
and charm quark coalescence will be important for J/ ip production, the J/ ip yield will be 
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Figure 5.9: J/ip yields from models as a function of number of participants. The black 

dashed curve depicts a calculation only includeing J/ip absorption effects, and no 
recombination effects from |Gra04j with cXcc = 0.35 mb. The black solid curve depict a 
calculation from statistical coalescence model and there is no absorption effect in 
it |ABR03j with acs = 0.39 mb. The dot-dashed lines depict several calculations 
including both absorption and recombination effects with different cross sections as 

input |Gra04j . 

"enhanced" in this case. Furthermore, the J/ip from coalescence production will have a 
larger {px) than those from direct production. Therefore, both enhancements in the J/tp 
yield and (pr) could be a good signature for QGP production at RHIC. The combination 
of the suppression and possible coalescence enhancement will, in principle, further enforce 
this argument. However, the energy loss of initially produced heavy quarks traversing 
the medium might complicate this picture. Detailed analysis of J/tp yields as well as 
transverse momentum distributions are necessary in order to distangle the problem. 

5.3 Elliptic flow of charm quarks and thermalization 

In heavy ion collisions, the momentum-space azimuthal anisotropy can reveal the early 
stage system information |ZGK99j . The large V2 measured for multi-strange hadrons and 
the NCQ scaling of hadrons at intermediate pt indicate a partonic level collectivity has 
been developed in Au -|- Au collisions at RHIC. However, partonic collectivity is necessary 
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but not sufficient to show that local thermal equilibrium has been established. Parton 
thermalization is crucial to demonstrate the discovery of QGP in heavy ion collisions at 
RHIC. 

Charm quarks are believed to be mostly produced through initial scatterings. Hence 
the study of elliptic flow of charm quarks may give us some hints of properties of the pro- 
duced matter in heavy ion collisions. Charm quarks, because of their heavy masses, need 
many more rescatterings to develop flow similar to light quarks {u, d, s). If charm quarks 
flow as much as light quarks, this may tell us there are frequent rescatterings happening 
between the light quarks and thus provide a strong clue to illustrate the thermalization 
of the light (u, d, s) quarks [DESoij. 
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Figure 5.10: V2 of electrons from the decays of charm hadrons and 7r°. Lines depict the 
input D or 7r° V2 while the points depict the V2 of decay electrons. Upper panel: the 
input D t>2 are assumed as V2 ~ ^2'^^ (dot-dashed line) and f 2 ~ (solid line), 
respectively. Bottom panel: the input f2(vr°) (dot-dashed line) is from 200 GeV 

measurement. 



At present, the direct reconstruction of open charm hadrons uses a mixed-event 
method. In this method, there are huge combinatorial backgrounds under the signal 
peaks, which leads to large uncertainties to extract V2 from this method. A new way to 
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measure the charm quark V2 is to measure the f 2 of their semi-leptonic decayed electrons. 
Fig. 15.101 upper panel shows the input open charm hadron f 2 and decayed electron f 2 
when charm quarks flow or not. About 50 M D^s were input into PYTHIA to obtain the 
decayed electrons. The shape of the px spectrum is a power law function which we 
discussed before. The results show that, at pt above 2 GeV/c, the electrons will carry 
almost all the saturated f 2 of open charm hadrons |DES04j and the electron f 2 can dis- 
tinguish between V2 ~ ^'2''^ and ~ 0. This demonstrates that measuring non-photonic 
electron V2 for p^ above 2 GeV/c is a feasible and efficient way to measure the charm 
quark jljKdOal lUEMj . 

Neutral pion decay is the dominant background source in the single electron analysis. 
The two photon processes are followed by conversion into electrons in the detector, and 
the rate of conversion depends on detector materials. These background electrons can be 
rejected by a topological method Joh02j . The remaining dominant background source is 
vr" Dalitz decay, which can only be subtracted statistically. The bottom panel of Fig. 15.101 
shows the electron V2 from 7r° Dalitz decays. The background electron flow leads to large 
uncertainties in the statistical subtraction. 

Electrons from heavy flavor decays begin to dominate the electron spectrum above 
Pt ~ 3 GeV/c. With the knowledge of the pion yield, D-meson yield, the pion V2, and the 
electron V2, it is possible to extract the D-meson V2. These measurements can be made 
by both the PHENIX and STAR collaborations at RHIC. Direct photon V2 can also be 
measured with this method. 
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CHAPTER 6 



Outlook 

6.1 Detector upgrade proposals 

STAR has proposed two important sub-detector upgrades: a full barrel Time-Of-Flight 
(TOF) detector |SlA()4j and a Heavy Flavor Tracker (HFT) |STA()5j . 

The proposed full barrel TOF detector would surround the outer edge of the TPC, 
and cover — 1 < r/ < 1 and ~ 27r in azimuth. By using the recently developed technology 
- MRPC, the TOF system can achieve the required timing resolution < 100 ps and the 
required particle detecting efficiency > 95%. This will significantly improve the STAR 
PID capability: t^/K can be separated to 1.8 GeV/c and p/meson can be separated up 
to ~ 3 GeV/c. The TOF detector allows us to greatly reduce the integrated luminosity 
needed for key measurements (collectivity, elliptic fiow etc.), to significantly extend the 
Pt reach of resonance measurements, and to study large and small scale correlations and 
fiuctuations etc.. Good performance of the prototype TOF detector (TOFr) in Run III 
and Run IV strongly offers us confidence about the full coverage barrel TOF detector. 

The Heavy Flavor Tracker detector will be placed between the beam pipe and the 
SVT detector. It will retain the (pr, coverage of the current TPC and have two 
layers of silicon pixel detectors. The 10 /im vertex resolution provided by the HFT will 
allow us to reconstruct open charm and bottom hadrons directly and make high precision 
measurements of px spectra, particle ratios, azimuthal anisotropy etc.. With the help 
of the SVT, it can be used to remove photon conversion electrons, which will make the 
measurement of the low mass e+e~ spectrum possible. 

At RHIC, several key measurements such as jet quenching, bulk collective motion, and 



87 



partonic collectivity etc. have provided strong hints for the discovery of QGP. However, 
at least two remaining critical points need to be demonstrated: whether the system is 
thermalized and whether the chiral symmetry is restored? These two sub-detector upgrade 
will help us to understand these questions and the two following measurements should be 
investigated. 



6.2 Open charm measurements 




I , 1 , \ , \ , 1 , u 

2 4 6 8 10 

Transverse momentum Pj (GeV/c) 

Figure 6.1: Raa of charm hadrons (thin lines) and decayed electrons (thick lines) for 
several assumed D spectra in Au + Au central collisions. The baseline used in p + p 
collisions is from STAR d + Au measurement. The thin dashed line depicts the Raa of 
D with the same behavior as charged hadrons. The thin dot-dashed line is for D with 
Tfo = 160 MeV and (Pt) = 0.4 c blast wave behavior in low pt and with the same Raa 
as charged hadron in high px ■ The thin solid line is for D with the same BW 
parameters as the dot-dashed line in low px , but the total yield of D is assumed to 
obey Nbin scaling. While in high pt , the Raa is taken from dead-cone energy loss 
calculation |DG04aj . Thick lines depict the electron distributions from the 
corresponding D distribution with the same line style. 

In the previous discussion, we proposed the measurement of charm azimuthal anisotropy 
through its semi-leptonic decay channel. However, the decay kinematics can smear the 
spectra differences. So it is hard to study charm meson radial flow or hydrodynamic 
behavior through its semi-leptonic decay channel. Fig. 16.11 shows Raa of charm hadrons 
and decayed electrons for several assumed D spectra in Au + Au central collisions. The 
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difference between the solid black line and the dashed black line for D is visible and it 
can be attributed to the smaller collectivity of charm hadrons due to heavy masses. But 
this difference is washed out by the decay kinematics and it is hard to tell the magnitude 
of charm hadron collectivity from decayed electron spectra. 

At high pt, heavy quark may lose less energy due to the "dead-cone" effect |DKOH 
IDG04aj . A precise charm-decayed electron spectrum measurement may give us some 
hints from the suppression factor. But, the contribution to electrons from bottom decays 
may become significant above 3 GeV/c and bottom quarks are expected to lose much 
less energy. As such, the electron spectrum suppression at high pt is hard to interpret in 
terms of charm hadrons and/or bottom hadrons separately. So the direct charm spectrum 
measurement is crucial to answer the hydrodynamic behavior and "dead-cone" effect of 
charm hadrons. Together with those measurements for light hadrons, we can understand 
the properties of the hot dense matter created in heavy ion collisions. 

Measured light hadron {u, d, s) yields and their ratios, from AGS to RHIC energies, 
have been well described in statistical models |Bra01j . Precise measurement on charm 
hadron yields and their ratios can help us understand the charm quark chemistry and 
the kinematical equilibration property of charm quarks in the medium. Meanwhile, since 
the charm quark yield is quite high at RHIC energy, coalescence models predict signifi- 
cantly different ratios of different charm hadrons |RS03j . Furthermore, the standard J/ip 
suppression scenario may be modified when thermally produced J/ip become significant. 
Precise measurements on these sensitive probes: /D^, /D^ and J/ip /D^ are needed 
to distinguish different pictures |STA05j . Table. 16.11 shows the events needed to observe 
3(T and D'^ signals with different detector configurations. With the help of TOP and 
HPT detectors, these precise measurements become feasible. 

The development of elliptic flow requires frequent rescatterings between components. 
Identified particle elliptic flow measurements show strong evidence that a partonic level 
collectivity has been established in Au + Au collisions at RHIC energy. However, it is 
still not sufficient to prove that the system is thermalized. Charm quarks are abundantly 
produced at RHIC energies. Due to their large masses and expected smaller hadronic 
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Table 6.1: Number of events needed to observe 3a and L>+ signals in Au + Au 200 
GeV collisions with different detectors configurations. 









TPC+SVT 


12.6 M 


500 M (X^ + K+) 


TPC+SVT+TOF 


2.6 M 


100 M 


TPC+SVT+TOF+HFT 


5 K (Minuit) 


700 K (Minuit + 7r+) 



rescattering cross sections, non-zero charm quark fiow would provide a strong indication 
that light quarks (w, s) are thermalized. Although an electron V2 measurement can 
provide an indirect effective method, the technical systematics need to be understood very 
well before we can draw conclusions. A precise direct measurement of the charm hadron 
fiow up to intermediate pt can be achieved with the upgraded detectors TOF and HFT. 
Simulations show the statistical errors can reach ~ 2 — 3% at pt ~ 5 GeV/c in one year 
of Au + Au run. 

6.3 Low mass e^e~ spectrum 

Electromagnetic probes have great advantages to investigate the early stage informa- 
tion of the system compared to hadronic probes because they don't have strong interac- 
tions with the final state system. The di-electron invariant mass M^+e- distribution has 
been proposed as a unique tool to study the properties of vector mesons (p, a;, 0, J/tl) etc) 
in-medium effect during heavy ion collisions. This is directly connected with the QGP 
signature — chiral symmetry restoration by measuring the mass position change of these 
vector mesons. Once the broken chiral symmetry is restored, {ipip) becomes and the 
quark mass coupling with the Higgs field decreases to 0. This is a robust proof of the 
discovery of QGP. 

This measurement has been carried out in low energy Pb + Au collisions at SPS. 
CERES/NA45 studied the e+e" invariant mass spectrum with projectile beams of 40A 
GeV and 158A GeV. They observed an excess in the di-electron spectrum compared to 
light hadron decays and attributed it to in-medium modifications of the p meson due to 
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the restoration of approximate chiral symmetry |Ada03aj . However, the data points from 
SPS have large errors both statistically and systematically, so we cannot draw a strong 
conclusion from them so far. In the phase diagram shown in Fig. II. 4t the SPS machine 
sits near the crossover boundary of the phase transition into QGP. While the Tch is higher 
and /ifi is lower in RHIC, this signal should be easier to observe. Since almost all other 
RHIC measurements indicate the formation of QGP, a precise measurement of the low 
mass e~^e~ spectrum at RHIC energy is likely to be interesting. 

This measurement needs both detector upgrade. By combining TOF and dE/dx in 
the TPC, electrons can be identified. This has been demonstrated in Fig. 13.11 and 13.191 
Background electrons from photon conversion and light hadron decays were studied in 
the simulations. Fig. 16.21 shows the simulated electron pair invariant mass distributions 
for those pairs from different sources in ~ 1 M PYTHIA + GEANT simulation events. 
The background level in the interesting uo^ p mass region is shown in Table 
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Figure 6.2: Me+e~ of electron pairs from PYTHIA + GEANT simulations for those 
from all sources (black solid circles), photon conversion only (red open circles) and vr^,?] 

Dalitz decays (blue triangles). 

According to the material properties of the detectors and also from simulations, photon 
conversions mostly happen in the SVT, especially the first layer of the SVT (out of three). 
So requiring charged particle hits in the SVT and HFT will significantly remove the photon 
conversion electrons. Fig. 16.31 shows the electron distribution from photon conversions 
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Table 6.2: Electron pair contribution from background 





UJ 





Total bkgd 


10-V(25 MeV/c") 


2 X 10-V(25 MeV/c^) 


7r° Dalitz 


5 X 10-7(25 MeV/c^) 


5 X 10-7(25 MeV/c^) 



with a 25-TPC-hit cut only (black), a 25-TPC-hit and a 2-SVT-hit cut (red) and a 25- 
TPC-hit, a 2-SVT-liit cut and a 2-HFT-hit cut (blue). With the help of SVT and HFT, 
~ 98% of conversion electrons can be removed. In this case, the 7r° will become the 
dominant background source for primary electrons. 
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Figure 6.3: Electrons from photon conversion spectrum for a 25-TPC-hit and a 
2-SVT-hit cut (red) and a 25-TPC-hit, a 2-SVT-hit cut and a 2-HFT-hit cut (blue). 
Only ~ 2% of total electrons from photon conversion are left with combined 

TPC+SVT+HFT cut. 



Furthermore, we can estimate the number of events needed to observe the vector 
meson signals under different background levels. From previous measurements and model 
simulations, we assume that dN'^ /dy ^ 1 in 200 GeV p + p collisions and dN'^ /dy ~ 300 
in Au + Au collisions; u/n ^ 0.15, ^/tt ^ 0.02, ^ 10 MeV/c^, ^ 8 MeV/c^ 
from simulations; The matching and electron PID efficiency with TOF+dE/dx method 
is ~ 80% and the matching efficiency of the combined TPC+SVT+HFT is ~ 60%. Then 
to observe 3cr signals for the u and 0, the number of events of Au + Au collisions needed 
are shown in Table. 16.31 
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Table 6.3: Number of events needed to observe 3(7 u and signals in Au + Au 200 
GeV collisions with different detectors configurations. 







4> 


TPC+TOF 


7 M 


2 M 


TPC+TOF+SVT+HFT 


800 K 


150 K 



The measurements of vector mesons through electromagnetic decay channels can shed 
light on whether chiral symmetry which is broken in QCD vacuum is restored in the QGP, 
and furthermore indicate the origin of mass in nature. 
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APPENDIX A 



QCD Lagrangian 

A.l Notations 

We use the natural units throughout this thesis: c = h = = 1- 

In the following part, we define the contravariant vectors in for space-time coordinates 
and 4-momentum vector in the Minkowski space. 

= a;\ a;^ x^) = (t, x) = (t, x, y, z) (A.l) 



= (p,p\p%Pl = {E,p) = (E,Pt,P,) = (E,p,,py,p,) (A.2) 



The space-time metric tensor g^^ is: 



/i ^ 

0-10 

0-10 

\o -ly 



(A.3) 



The covariant vector is related to the contravariant vector through g^^, by 



(A.4) 



and conversely, 



(A.5) 
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In this thesis, we use the notation that a repeated index imphes a summation with respect 
to that index, unless indicated otherwise. The scaler product of two vectors a and h is 
defined as: 

a-h = a^h^ = g^^a'^h" = a°6° -a b (A.6) 
The gradient operator is the derivative with respect to x^: 

The four-momentum operator in coordinate representation is 

pM ^ idmu ^ .g,uQ^ ^ ^g,ud ^ d _^d _^d _^d 

ox^ ox" OX^ OX'^ OX"^ 

It is convenient to work with gamma matrices later. They satisfy the anticommutation 
relation 

{7^7'^} = 7V + 7V = 2/'^ (A.9) 

A. 2 SU{3)c invariant QCD Lagrangian 

We denote qj a quark field with color a and fiavor / and adopt a vector notation in 
color space: qf = colu'mn{qj, qj, q^). The free Lagrangian is written as: 

Co = J2qfitrd,-mf)qf (A.IO) 
/ 

It is invariant under arbitrary global SU{3)c transformations in color space: 

qf {q])' = [/^g^, UU^ = U^U = 1, detC/ = 1 (A.ll) 

The U (?>)c matrices can be written as 

U = exp (-ig.^e^ (A.12) 
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where A" (a — 1, 2, 8) denote the generators of the fundamental representation of the 
SU{?>)c algebra, and 9a are arbitrary parameters. In SU{?>)c, the matrices A" correspond 
to the eight Gell-Mann matrices and they satisfy the commutation relations 



[A", A"] = A'' A" - A' A" = 2ir''\^ (A. 13) 

with /"^'^ structure constants, which are real and totally antisymmetric. 

As in the QED case, to satisfy the invariance under local SU{3)c transformations: 
9a = 9a{x), 8 independent gauge bosons G'^{x) - gluons are introduced. Define: 



qf = [d^^-ig,G^^{x)]qf (A.14) 



To require D^^qf to transform in exactly the same way as g/, the transformation properties 
of the gauge fields under an infinitesimal SU{2))c transformation: 



G'a ^ {G^a)' = G^a- d^i^Oa) + 9sn5d,G^, (A.15) 

Unlike the QED case, the non-commutativity of the SU{2))c matrices gives rise to an 
additional term involving the gluon fields themselves. To build a gauge invariant kinetic 
term for the gluon fields, the field strengths were introduced: 



9s 2 



(A.16) 



The final SU{?>)c invariant QCD Lagrangian is: 

^QCD ^ -WgI^ + J2 <lAn'D, - ^fhf (A.17) 
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It is worthwhile to decompose the Lagrangian into different pieces: 

f 

f \ ^ / a/3 

_ |/afc(a^G^ - d'^G^jG'^Gl - 3lfa^cj^Q^^Q.Qpe 

The first hne contains the correct kinetic terms for different fields. The second line shows 
the color interaction between quarks and gluons. In the last line, owing to the non-abelian 
character of the SU{Z)c group, the G'^'^G'^^^^ term generates the cubic and quartic gluon 
self- interactions. 

A. 3 Chiral symmetry and effective chiral Lagrangian 

In the absence of quark masses, the QCD Lagrangian can be written: 

Cqcd = -\g7G% + mrD^QL + iqRrD.QR (A.19) 

. It is invariant under independent global G = SU{Nf)L <8) SU{Nf)R transformations of 
the left-and right-handed quarks in flavor space: 

G G 

Ql gLQL, Qr ^ QrQr, 9l,r e SU{Nf)L^R (A.20) 

This chiral symmetry, which should be approximately good for light quark sector 
{u, d, s), is however not seen in the hadronic spectrum. Moreover, the octet of pseudoscalar 
mesons are much lighter than all the other hadronic states. Hence, the ground state 
of the theory (the vacuum) should not be symmetric under the chiral group. In the 
Goldstone's theorem, an octet of pseudoscalar massless bosons is introduced and the 
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symmetry SU{3)l ^ SU{3)r spontaneously breaks down and the quark condensate 



V = {0\uu\0) = {0\dd\0) = (0|ss|0) < 



(A.21) 



The Goldstone nature of the pseudoscalar mesons imphes strong constraints on their 
interactions, which can be most easily analyzed on the basis of an effective Lagrangian. 
The Goldstone boson fields can be collected in a 3 x 3 unitary matrix U{(f)). 



(A.22) 



A convenient parametrization is given by 



U{(f)) = exp(iV2$//) 



(A.23) 



~2 



— + ^ 

V2 



V 



TT 



K- 



K+ 



\ 



V6j 



(A.24) 



The matrix U{(f)) transforms linearly under the chiral group, but the induced transforma- 
tion on the Goldstone fields (p is highly non-linear. 

We should write the most general Lagrangian involving the matrix U (0) and organize 
it in terms of increasing powers of the momentum or, equivalently, in terms of an increasing 
number of derivatives: 

Ceff{U)=J2^2n (A.25) 



The terms with a minimum number of derivatives will dominate in the low energy domain. 
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To lowest order, the effective chiral Lagrangian is uniquely given by the term 




(A.26) 



The non-linearity of the effective Lagrangian relates amplitudes with different numbers 
of Goldstone bosons, allowing for absolute predictions in terms of /. 

Considering the quark masses, the corrections induced by the non-zero masses are 
taken into account through the term 



The relation between the physical meson masses and the quark masses can be obtained 
rom the trace in the quadratic mass term. Although the absolute values of the quark 
masses cannot be fixed from this approach because of the factor \v\/ f in each mass, one 
can obtain information about quark mass ratios after taking out this common factor. One 
of the famous ratios advocated by Weinberg is: 




(A.27) 



m„: m^: rUg — 0.55: 1: 20.3 



(A.28) 
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APPENDIX B 
Kinematic variables 



Let us introduce several useful kinematic variables in hadron-hadron interactions. 

Pz (^) 

• Bjorken x, Xsjorken = — 7T~~, 7- the longitudinal momentum fraction carried by 

Pz{hadron) 

the parton i to the total hadron momentum. Usually denoted as x. 

p* 

• Feynman x, xp = — — r: the longitudinal momentum fraction of final particle 

p*[max) 

to the maximum momentum in the center-of-mass system. 

• transverse mass rriT = y^pfT+m^' 

• rapidity y = - In [ ^° ] . Then po — uit cosh y, and Pz — rriT sinh y. 



2 Vi^o - Pz 

• pseudorapidity r? = - In ( l^l— — ^ ) • Then Ipl = pt cosh, ri, and ^p^sinhTi. 

2 \\P\-Pzj 

We consider a two-parton interaction 1 + 2 ^ 1' + 2' + .... Two partons are from the 
incoming hadron beams of symmetric energy ^/s . Usually, the z axis is defined to parallel 
to the beam. At high energy {^Js ^ mh) , before interaction, the four- momentum of the 
two partons are: 

p'l = {x,^, 0, 0, x,^), pi^ = 0, 0, -X2^) (B.l) 

The four-momentum of final state particle j is: 

pf-ip'opPT'pP'.j) (B.2) 
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Due the momentum conservation 



' 'I'l 'I'n I > 7^ 



Using p'q- = p'rp- cosh|/j and p' • = p^ - siniiyj, we can obtain 



1 



3 



(B.4) 



If we consider a two-particle final state, transverse momentum conservation requires = 
Pt2 — Pt (neglect the symbol /), at mid-rapidity where yi^2 ~ 0, we have: 

xi = X2 = xt ^ — ^ (B.5) 

In general cases, we often use this treatment XBjorken ^ xt — 2pT/y/s. 

Since the maximum fraction of the longitudinal momentum for a certain parton is 
the Feynman x is easily to be extracted: 

XF^^ (B.6) 
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APPENDIX C 



Low energy charm cross section data points selection 

There were a lot of charm cross section measurements at low energy {^/s < 70 GeV). 
Table O lists all the charm cross section data points at low energies from journal publi- 
cations. We ignored those data points extrapolated from high xp and/or extremely low 
efficiency from correlation measurements, with extrapolation factor of ^ 10. The refer- 
ences for these data points are listed as followings. The data values included in the thesis 
are multiplied by 1.5 to account for the additional Df and contribution |FMN97j . 



1. S.P.K. Tavernier, Rep. Prog. Phys. 50, 1439(1987) and references therein. 

2. NA32 collaboration, Z. Phys. C39, 451(1988). 

3. E769 collaboration, Phys. Rev. Lett. 77, 2388(1996). 

4. NA16 collaboration, Z. Phys. C40, 321(1988). 

5. NA27 collaboration, Phys. Lett. B 135, 237(1984). 

6. E743 collaboration, Phys. Rev. Lett. 61, 2185(1988). 

7. E653 collaboration, Phys. Lett. B 263, 573(1991). 
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Table C.l: Overview of charm cross section data points at low energies and the criteria 
to include those in the thesis. Some experiments used da/dxp = (1 — for 

extrapolation. 









V^'UllilliClll 


Included? 


(GeV/c) 


+ Target 






Ref. 


c\r\r\ 

200 


1 o * 

p + bi 


1.5 ± 0.7 


NA32,iJ", iJ^,good 
acceptance (0-20%) at > 


/ 


[2] 


200 


p + CsFs 


3-9Ii.9 


NA25, Bubble chamber HOBC, 1// 
trigger, > o CeV, good 
acceptance at Xi? > (8% efficiency) 


/ 


[1] 


250 


p + Be,Uu, 
A1,W 


y.o ± 1.5 


i!j769, D\ U^, good acceptance 
at xf > 


/ 


[3] 


2(50 


„ 1 TT^ /IVT^ 

n + hLe/lNe 


5.6 ± 1.7 


E630 streamer chamber, 
1/U trigger, £'(/x) >6.5 GcV 


/ 


[1] 


ooO 


p -\- be 


22 ± y 


bll-btaniorcl, prompt Ifi trigger, 
-C/(yUj >zU LreV, acceptance 09/0 


/ 


[1] 


OOU 


= 

p -\- re 


11 Q _L 9 n 


r iN/\ij-v^v^r rv,o, pionipt iyu. tiiggei, 
E(u) >20 CeV Tl— 5 r >0 3 


/ 


[1] 


360 


P + CgFg 


24.61^:^3^ 


See 200GeV data 


/ 


[1] 


360 


pn- p 


'J-'J-4.6 


NA16 Ho bubble chambpr T,FRC-I-F,HS 
good acceptance for xj? >0 


/ 


r ^ 1 

[4] 


400 


p -t- p 


1 c, 1 _l_ 1 7 


NA97 H9 bnbblp rhpmbpr 

good acceptance for >0 




r r 1 

[5] 


400 




"5 -"--18 


CTT-Stanfnrd 1;/ F,( ii) ->90 QpV 
1.0 < Pt < 2.5 GeV, 2.5%, n~5 


/ 


[1] 


soo 


P -t- P 


8 


F743 T.FRC-MPS -Ol^rr^^O^ n~8 




[6] 


800 


P + P 


'"-19 


E653, -0.2< xp <0.5, n~7 




[7] 


/i c\c\ 

400 


p + Cu 




GGM 1/ beam dump, n=3, [x) =0.8 


X 


[1] 


400 


p + Cu 


on 1 n 
39 ± 10 


CFiKiN 1st u beam dump, CDBb, 
■ty\fJ>) > .^u vjev, n — o-o, \x) — u.o 


X 


[1] 


400 


p + Cu 


17 ±4 


CERN-BEBC, 2nd beam dump, 
E(/i) > 10 GcV, n=4, (x) =0.8 


X 


[1] 


400 


p + Cu 


15 ±5 


CERN-CHARM u beam dump, 
E{ij) > 20 GeV, n=4, {x) =0.8 


X 


[1] 


400 


p + W 


15.5 ±2.5 


E613 FNAL u beam dump, 
E{iJ,) > 20 GeV, n=4, (x) =0.45 


X 


[1] 


400 


p + Cu 


~ 17 


WA66-BEBC, CERN 2nd u beam dump, 
E(^) > 10 GcV, n~5 


X 


[1] 


400 


p + Cu 


15.5 ±2.9 


CHARM, CERN 2nd ly beam dump, n=5 


X 


[1] 


400 


p + Cu 


~20 


CDHS, CERN 2nd u beam dump, n=5 


X 


[1] 



the aboA'c experiiiieiits co\'er high x. large extrapolation errors 
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continued 



Beam p 


Beam 


(^cc (/wb) 


Comment 


Included? 


(GeV/c) 


-|- Target 






Ref. 


400 


p -|- Fe 


7~ 20 


CIT-Stanford, 2/i-|-missing energy, 
Pt >0.75 GeV, 0.1-0.4% acceptance 


X 


[1] 


400 


p + W 


< 10 


H + H mass 


X 


[1] 


using correlation, low acceptance 












63 


p + p 


840 ± 320 


ISR, trigger, Kc{K-np) only, n~0 


X 


[1] 


63 


p + p 


150 - 450 


ISR, e trigger, kciKup) and D{Knn) 


X 


[1] 


62 


p + p 


650 ± 222 


ISR, CBF, e trigger, n=3(L'), n=0(Ae) 


X 


[1] 


62 


p + p 


129 ± 75 


ISR, CBF, new runs, Kc{K'Kp) only 


X 


[1] 


53-62 


p + p 


1390 ± 180 


ISR, Kc{KTip) only, 
n=0, X >0.75 


X 


[1] 




correlation, hij 


Xpi large extrapolation errors 






53-62 


p + p 


70 ±36 


ISR, e + 11 unlike-sign pairs, 
0.0016% efficiency, n~5 


X 


[1] 


53-62 


p + p 


73 ± 21 


ISR. cc pairs 


X 


[1] 



correlation, low efficiency 
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APPENDIX D 



Presentations and publication list 
Presentations 

• Open charm production at RHIC - recent results from STAR 

21st Winter Workshop on Nuclear Dynamics, Breckenridge, Colorado, USA, 02/05/2005 
- 02/12/2005. 

• Elliptic flow of pion, kaon, proton from Au + An collisions at 62.4 GeV 

2004 Fall Meeting of the Division of Nuclear Physics of APS, Chicago, Illinois, USA, 
10/27/2004 - 10/31/2004. 

• Open Charm Yields in d + Au Collisions at 200 GeF (poster) 

2004 Gordon Research Conference on Nuclear Chemistry, New London, New Hamp- 
shire, USA, 06/13/2004 - 06/18/2004. 

• Resonance decay effects on Anisotropy Parameters 

April Meeting of APS, Denver, Colorado, USA, 05/01/2004 - 05/04/2004. 

• Open Charm Yield in 200 GeV d + Au Collisions at RHIC 

2004 April Meeting of APS, Denver, Colorado, USA, 05/01/2004 - 05/04/2004. 

• The Performance of a Prototype Multigap Resistive Plate Chamber Time- Of- Flight 
Detector for the STAR Experiment (poster) 

Quark Matter 2004, Oakland, California, USA, 01/11/2004 - 01/17/2004. 

• Single Electron Spectra from d-hAu and p-hp collisions at y^sjviv — 200 GeV 
2003 Fall Meeting of the Division of Nuclear Physics of APS, Tucson, Arizona, USA, 
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10/28/2003 - 10/31/2003. 



Publication List 

• Open charm yields in d + Au collisions at -s/snn — 200 Ge V 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 94, 062301(2005) 
Principle authors: X. Dong, L. Ruan, Z. Xu and H. Zhang. 

• Resonance decay effects on anisotropy parameters 

X. Dong, S. Esumi, R Sorensen, N. Xu and Z. Xu, Phys. Lett. B597, 328(2004). 

• Open charm production at RHIC - recent results from STAR 

X. Dong et al. (for STAR Collaboration), Proceedings of 21st Winter Workshop on 
Nuclear Dynamics. 

• Improvement on the charge resolution with the average of truncated dE/dx and op- 
timized combination 

X. Dong, S.W. Ye, H.F. Chen, Z.P. Zhang and Z.Z. Xu, Journal of University of 
Science and Technology of China Vol 8 (2002) (in Chinese). 

• Pion, kaon, proton and anti-proton transverse momentum distributions from p -\-p 
and d + Au collisions at ^s^n = 200 GeV 

J. Adams et al. (STAR Collaboration), arXiv: nucl-ex/ 0309 102 
Principle authors: L. Ruan, X. Dong, F. Geurts, J. Wu and Z. Xu. 

• Azimuthal anisotropy and correlations at large transverse momenta inp-\-p and Au 
-h Au collisions at -y/sNN — 200 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 93, 252301(2004). 

• Azimuthally sensitive HBT in Au -h Au collisions at ^sjvat — 200 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 93, 012301(2004). 

• Multi-strange baryon production in Au-Au collisions at sJsnn — 130 GeV 
J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 182301(2004). 
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Cross Sections and Transverse Single-Spin Asymmetries in Forward Neutral Pion 
Production from Proton Collisions at y/s — 200 Ge V 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 171801(2004). 

Identified particle distributions in pp and Au + Au collisions at a/s^n — 200 GeV 
J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 112301(2004). 

p° Production and Possible Modification in Au + Au andp + p Collisions at y^Sivjv 
= 200 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 092301(2004). 

Azimuthal anisotropy at the Relativistic Heavy Ion Collider: the first and fourth 

harmonics 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 062301(2004). 

Particle-type dependence of azimuthal anisotropy and nuclear modification of particle 
production in Au-hAu collisions at sJsnn — 200 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 052302(2004). 

Pion-Kaon Correlations in Central Au -h Au Collisions at sJsnn — 130 GeV 
J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 91, 262302(2003). 

Three-Pion Hanbury Brown- Twiss Correlations in Relativistic Heavy-Ion Collisions 
from the STAR Experiment 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 91, 262301(2003). 

Transverse momentum and collision energy dependence of high px hadron suppres- 
sion in Au+Au collisions at ultrarelativistic energies 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 91, 172302(2003). 

Evidence from d+Au measurements for final-state suppression of high px hadrons 
in Au -h Au collisions at RHIC 

J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 91, 072304(2003). 
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Pseudorapidity Asymmetry and Centrality Dependence of Charged Hadron Spectra 
in d + Au Collisions at sJsnn) — 200 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 064907(2004). 

Measurements of transverse energy distributions in Au + Au collisions at s/smn — 
200 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 054907(2004). 

Photon and neutral pion production in Au + Au collisions at s/snn — 130 GeV 
J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 044902(2004). 

Centrality and pseudorapidity dependence of charged hadron production at interme- 
diate pt in Au + Au collisions at ^/sNN = 130 GeV 
J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 044901(2004). 

Rapidity and Centrality Dependence of Proton and Anti-proton Production from Au 
-h Au Collisions at ■•^snn — 130 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 041901(2004). 

Production of e'^e~ Pairs Accompanied by Nuclear Dissociation in Ultra- Peripheral 
Heavy Ion Collision 

J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 031902 (R) (2004). 

Transverse-momentum dependent modification of dynamic texture in central Au -h 
Au collisions at ^s^^ = 200 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. C 70, 031901 (R) (2004). 

Net charge fluctuations in Au + Au collisions at y^sjvTV = 130 GeV 

J. Adams et al. (STAR Collaboration), Phys. Rev. C 68, 044905(2003). 
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